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1 Introduction 



Topological defects are ubiquitous in physics. Whenever a sym- 
metry breaking phase transition occurs, topological defects may 
form. The best known examples are vortex lines in type II super 
conductors or in liquid Helium, and declination lines in liquid 
crystals [ |iU9| , |2^] . In an adiabatically expanding universe which 
cools down from a very hot initial state, it is quite natural to pos- 
tulate that topological defects may have emerged during a phase 
transition in the early universe and that they may have played 
the role of initial inhomogeneities seeding the formation of cosmic 
structure. This basic idea goes back to Kibble (1976) In this 
report we summarize the progress made in the investigation of 
Kibble's idea during the last 25 years. Our understanding of the 
formation and evolution of topological defects is reported almost 
completely in the beautiful book by Vilenkin & Shellard \\ib'Z 



or 



the excellent Review by Hindmarsh & Kibble [j7T|, and we shall 



hence be rather short on that topic. Nevertheless, in order to be 
self contained, we have included a short chapter on spontaneous 
symmetry breaking and defect formation. Our main topic is how- 
ever the calculation of structure formation with defects, results 
which are not included in [|152|j and [)7T[ . 



Besides the formation of structure in the universe, topological 
defects may be relevant for the baryon number asymmetry of the 
|33|j . Superconducting cosmic strings [p^j or vortons [pO| 



universe 



might produce the high energy cosmic rays [jT5j], or even gamma 
ray bursts [jT2j]. The brane worlds which have have focussed a 
lot of attention recently, may actually just represent topological 
defects in a higher dimensional space [|55| , |6T| , |62|] . There have also 
been interesting results on chiral strings and their observational 
signatures [PT| , |i38|] . GUT scale cosmic strings could be detected 
by their very peculiar lensing properties. For a straight cosmic 
string lensing is very simple [ |152| . For a more realistic network 
of strings, characteristic caustics and cusps in the lensing signal 
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are very generically expected [p6| , |i4«| . 



The relevant energy scale for a topological defect is the phase 
transition temperature. Hence a good estimate for the amplitude 
of the dimensionless gravitational potential induced by topo- 
logical defects is 

* - 47rGT^ = 47r{TjMpf , (1) 



where Mp denotes the Planck mass. The measurements of cos- 
mic microwave background anisotropics on large scales by the 
cosmic background explorer (COBE) satellite [|ia4j] have found 



that this potential, which is of the same order as the tempera- 
ture fluctuations on large scales, is about 10"^. Hence, for cos- 
mic structure formation, we are interested in phase transitions at 
Tc ~ 10~^Mp ~ lO^^GeV. Interestingly, this is just the scale of 
the grand unification phase transition (GUT scale) of supersym- 
metric GUT's (grand unified theories). 

Topological defects represent regions in space-time where the cor- 
responding field (order parameter in condensed matter physics or 
Higgs field in particle physics) is frustrated. It cannot relax into 
the vacuum state, the lowest energy state, by topological obstruc- 
tions. They represent positions of higher energy density and are 
thus inherently inhomogeneous distributions of energy and mo- 
mentum. We shall discuss concrete examples later. 

In the remainder of this introduction we give a brief overview 
of the problem of structure formation and we present the main 
results worked out in this report. 

In Chapter g we introduce the concept of topological defect for- 
mation during symmetry breaking phase transitions, we classify 
the defects and illustrate them with examples. 
In Chapter we present in detail the theoretical framework used 
to investigate structure formation with topological defects. This 
chapter together with two appendices is self contained and should 
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enable a non-specialist in the field to fully understand the often 
somewhat sketchy literature. 

In Chapter g we discuss numerical simulation of topological de- 
fects. We distinguish global and local defects which have to be 
treated in a very different way. We specify the approximations 
made in different numerical implementations and discuss their 
validity and drawbacks. 

In Chapter |^ we present the results of simulations of structure for- 
mation with topological defects and compare them with present 
observations. 

In Chapter ^ we investigate the question in whether the results 
discussed in Chapter |^ are generic or whether they are just spe- 
cial cases. We derive properties of the unequal time correlators of 
generic causal scaling seeds. Since these are the sole ingredients in 
the calculation of the fluctuation power spectra, they determine 
the 'phase space' of defect models of structure formation. We 
discuss a model of causal scaling seeds which mimics the cosmic 
microwave background (CMB) anisotropy spectrum of inflation. 
We also consider the possibihty that large scale structure may be 
due to a mixture of adiabatic inflationary initial perturbations 
and topological defects. We study especially the fluctuations in 
the CMB. We investigate to which extent CMB parameter esti- 
mations are degraded if we allow for an admixture of defects. 

We end with a brief summary of the main results. 

Throughout this work we use units with c = h = /cBoitzmann = 1- 
The basic unit is thus either an energy (we usually take MeV's) 
or a length, for which we take cm or Mpc depending on the sit- 
uation. 

We choose the metric signature (— , +, +, +). Three-dimensional 
vectors are denoted in boldface. The variables x and k are comov- 
ing position and comoving wave vector in Fourier space. Greek 
indices, /i, z/, • • • denote spacetime components of vectors and ten- 
sors while Latin indices j, • • • denote three dimensional spatial 
components. We mostly use conformal time r with dr = dt/a, 
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where t is cosmic time and a is the scale factor. Derivatives with 
respect to conformal time are indicated by an over-dot, f = ^■ 



1.1 Main results 



Before we start to discuss models of structure formation with 
topological defects in any detail, let us present the main results 
discussed in this review. 

We concentrate primarily on CMB anisotropics. Since these ani- 
sotropics are small, they can be calculated (almost fully) within 
linear cosmological perturbation theory. To compare models with 
other data of cosmic structure, like the galaxy distribution or ve- 
locities, one has to make some assumptions concerning the not 
well understood relation between the distribution of galaxies and 
of the dark matter, the problem of biasing. Furthermore, on small 
scales one has to study non-linear Newtonian clustering which is 
usually done by A^^-body simulations. But to lay down the initial 
conditions for A^-body simulations, one does not only need to 
know the linear power spectrum, but also the statistical distribu- 
tion of the fluctuations which is largely unknown for topological 
defects. Fluctuations induced by topological defects are generi- 
cally non-Gaussian, but to which extent and how to characterize 
their non-Gaussianity is still an open question. In this report, 
we therefore concentrate on CMB anisotropics and their polar- 
ization and shall only mention on the side the induced matter 
fluctuations and bulk velocities on large scales. 

Like inflationary perturbations, topological defects predict a Har- 
rison-Zel'dovich spectrum of perturbations [|6^ , |it)a|] . Therefore, 



the fluctuation spectrum is in good agreement with the COBE 
DMR experiment [ |i34|] , which has measured CMB anisotropics 



on large angular scales and found that they are approximately 
constant as expected from a Harrison-Zel'dovich spectrum of ini- 
tial fluctuations (see e.g. [|TT3]). 



7 



Since quite some time it is known, however, that topological 
defect differ from adiabatic inflationary models in the acoustic 
peaks of the CMB anisotropy spectrum [Q. Due to the isocur- 



vature nature of defects, the position of the first acoustic peak is 
shifted from an angular harmonic of about £ ~ 220 to £ ~ 350 
(or up to £ ~ 500 in the case of cosmic strings) for a spatially 
flat universe. More important, the peaks are much lower in de- 
fect models and they are smeared out into one broad hump with 
small wiggles at best. Even by changing cosmological parameters 
at will, this second characteristics cannot be brought in agree- 
ment with present CMB anisotropy data hke [|ii0| , |55| , |g7|] . Also 



the large scales bulk velocities, which measure fluctuation ampli- 
tudes on similar scales turn out to be too small [|57| . 



As the CMB anisotropy signals from cosmic strings and from 
global 0{N) defects are quite different, it is natural to wonder 
how generic these results may be. Interestingly enough, as we 
shall see in Chapter g, one can define 'scaling causal seeds', i.e. 
initial perturbations which obey the basic constraints for topolog- 
ical defects and which show a CMB anisotropy spectrum resem- 
bling adiabatic inflationary predictions very closely [|i44|] . This 



'Turok model' can nevertheless be distinguished from adiabatic 
perturbations by the CMB polarization spectrum. Also mixed 
models with a relatively high degree of defect admixture, up to 
more than 50%, are in good agreement with the data. 



1.2 Cosmic structure formation 



The geometry of our universe is to a very good approximation 
isotropic and therefore (if we assume that we are not situated 
in a special position) also homogeneous. The best observational 
evidence for this fact is the isotropy of the cosmic microwave 
background which is (apart from the dipole anisotropy) on the 
level of about 10"^ - 10^1 

Nevertheless, on galaxy and cluster scales, the matter distribu- 
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tion is very inhomogeneous. If these inhomogeneities have grown 
by gravitational instabihty from small initial fluctuations, the 
amplitude of the initial density fluctuations have to be about 
10^^ to become of order unity or larger today. Radiation pres- 
sure inhibits growth of perturbations as long as the universe is 
radiation dominated, and even in a matter dominated universe 
small perturbations grow only like the scale factor. 

The discovery of large angular scale fluctuations in the CMB by 
the DMR experiment aboard the COBE satellite [ |ia4|] , which are 
just about of the needed amplitude, is an important support of 
the gravitational instability picture. The DMR experiment also 
revealed that the spectrum of fluctuations is 'flat', which means 
that the fluctuations have a fixed amplitude A when entering the 
Hubble horizon. This implies that temperature fluctuations on 
large scales are constant, as measured by COBE: 

^(^^j (i?)^ ~ IQ-if^ = constant, (2) 

independent of the angle i!}, for i!) ^ 1°. These fluctuations are 
of the same order of magnitude as the gravitational potential. 
Their smallness therefore justifies the use of linear perturbation 
theory. Within linear perturbation theory, the gravitational po- 
tential does not grow. This observation originally led Lifshitz to 
abandon gravitational instability as the origin for cosmic struc- 
ture [I'iUUlj . But density fluctuations of dust can become large and 
can collapse to form galaxies and even black holes. At late times 
and on scales which are small compared to the horizon, linear per- 
turbation theory is no longer valid and numerical N-body simu- 
lations have to be performed to follow the formation of structure. 
But also with N-body simulations one cannot compute the details 
of galaxy formation which strongly depend on thermal processes 
like cooling and formation of heavy elements in stars. Therefore, 
the relation of the power spectrum obtained by N-body simula- 
tions to the observed galaxy power spectrum may not be straight 
forward. This problem is known under the name of 'biasing'. 
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Within linear perturbation analysis, structure formation is de- 
scribed by an equation of the form 



where "D is a time dependent linear differential operator, k is 
the wave vector and X is a long vector describing all the cosmic 
perturbation variables for a given fc-mode, like the dark matter 
density and velocity fluctuations, the CMB anisotropy multipole 
amplitudes and so on. <S is a source term. In Chapter ^ we will 
write down the system (|^) explicitly. 

There are two basically different classes of models for structure 
formation. In the first class, the linear perturbation equations 
are homogeneous, i.e. S = 0, and the resulting structure is de- 
termined by the initial conditions and by the cosmological pa- 
rameters of the background universe alone. Inflationary initial 
perturbations are of this type. In most models of inflation, the 
initial fluctuations are even Gaussian and the entire statistical 
information is contained in the initial power spectrum. The evo- 
lution is given by the differential operator T> which depends on 
the cosmological parameters. 

In the second class, the linear perturbation equations are inho- 
mogeneous, having a so called source term or 'seed', <S, on the 
right hand side which evolves non-linearly in time. Topological 
defects are of this second class. The difficulty of such models hes 
in the non-linear evolution of the seeds, which in most cases has 
to be determined by numerical simulations. Without additional 
symmetries, like e.g. scaling in the case of topological defects, it 
is in general not possible to simulate the seed evolution over the 
long time-scale and for the considerable dynamical range needed 
in cosmology. We shall see in Chapter g how this problem is over- 
come in the case of topological defects. An additional difficulty 
of models with seeds is their non-Gaussian nature. Due to non- 
linear evolution, even if the initial conditions are Gaussian, the 
seeds are in general not Gaussian at late times. The fluctuation 




(3) 
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power spectra therefore do not contain the full information. All 
the reduced higher moments can be relevant. Unfortunately only 
very little work on these non-Gaussian aspects of seed models has 
been published and the results are highly incomplete [|6D| , ^ , ^ 



11 



2 Symmetry Breaking Phase Transitions and the Formation of 
Topological Defects 

2.1 Spontaneous symmetry breaking 

Spontaneous symmetry breaking is a concept which originated in 
condensed matter physics. As an example consider the isotropic 
model of a ferro-magnet: although the Hamiltonian is rotationally 
invariant, the ground state is not. The magnetic moments point 
all in the same direction. 

In models of elementary particle physics, symmetry breaking is 
most often described in terms of a scalar field, the Higgs field. In 
condensed matter physics this field is called the order parameter. 
It can also be a vector or tensor field. 

A symmetry is called spontaneously broken if the ground state 
is not invariant under the full symmetry of the Lagrangian (or 
Hamiltonian) density. Since the symmetry group can be repre- 
sented as a group of linear transformations, this implies that the 
vacuum expectation value of the Higgs field is non-zero. 

The essential features of a spontaneously broken symmetry can 
be illustrated with a simple model which was first studied by 
Goldstone (1961) [Q. This model has the classical Lagrangian 



Here is a complex scalar field and A and t] are real positive con- 
stants. The potential in (g) is called the 'Mexican hat potential' 
as it looks somewhat like a Mexican sombrero. This Lagrangian 
density is invariant under the group U{1) of global phase trans- 
formations. 



density 






(5) 
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The minima of the potential V lie on the circle |0| = r/ which is 
called the 'vacuum manifold', A4 = S,^^ (here and in what follows 
denotes an n-sphere of radius R and S" denotes an n-sphere 
of radius 1). The notion 'global' indicates that the symmetry 
transformation is global, i.e., a is independent of the spacetime 
position X. The quantum ground states (vacuum states) |0) of 
the model are characterized by 

(O|0|O) = r]e"^ ^0 . (6) 



A phase transformation changes (3 into + a, hence a ground 
state is not invariant under the symmetry transformation given 
in Eq. (|^). (Clearly, the full vacuum manifold A4 is invariant 
under symmetry transformations and thus a mixed state which 
represents a homogeneous mixture of all vacuum states is still 
invariant even though no pure state is.) The only state \u) in- 
variant under the symmetry (|5]), characterized by = 0, 
corresponds to a local maximum of the potential. Small pertur- 
bations around this 'false vacuum' have 'negative mass' which 
indicates the instability of this state: 

V(0) = -^Ar/2|0|2 + const. + O(|0n . (7) 

The vacuum states of the broken symmetry are all equivalent and 
we can thus reveal their properties by studying one of them. For 
convenience we discuss the vacuum state with vanishing phase, 
(O|0|O) = T]. Expanding the field around this state yields 

0(x) = (r/ + -i=^(x))e^^(-), (8) 

where (p and 'd are real fields. The Lagrangian density in terms 
of (p and 'd becomes 

j0=1 {d,^f + rf {d,df - \\T^V + CUv, . (9) 
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The interaction Lagrangian Cmt is easily determined from the 
original Lagrangian, (|^). This form of the Lagrangian shows that 
the degree of freedom if is massive with mass = Xif while 'd 
describes a massless particle (circular excitations), a Goldstone 
boson. This simple model is very generic: whenever a continu- 
ous global symmetry is spontaneously broken, massless Goldstone 
bosons emerge. Their number equals the dimension of the vacuum 
manifold, ie., the dimension of a group orbit (in the space of field 
values). In our case the space of field values is C « R^. A group 
orbit is a circle of dimension 1 leading to one massless boson, the 
excitations tangential to the circle which cost no potential energy. 
The general result can be formulated in the following theorem: 

Theorem 1 (Goldstone, 1961) HU^] If a continuous global sym- 
metry, described by a symmetry group G is spontaneously broken 
to a sub-group H C G, massless particles emerge. Their num- 
ber is equal to the dimension n of the vacuum manifold A4 ( the 
"number of broken symmetries"). Generically, 

M. = G/ H and n = diniG — dimH = dim {G/H) , 

where here = means topological equivalence. 

In our example G = U{1), H = {1} and n = 1 — = 1. 
Another well-known example are the three pions, tt*, tt*^, which 
are the Goldstone bosons of isospin (proton/neutron) symmetry. 
There the original symmetry, SU (2) is completely broken leading 
to n = dimSU{2) = 3 Goldstone bosons (see. e.g. [[TTl ). 

Very often, symmetries in particle physics are gauged (local). The 
simplest gauge theory is the Abehan Higgs model (sometimes also 
called scalar electrodynamics). It is described by the Lagrangian 
density 

jC = D.^D'^c^ - y(0) - \f,,F^'' , (10) 
where is again a complex scalar field and Dn = dn — ieAn is the 
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covariant derivative w.r.t. the gauge field A^^. = df^A,^ — d„Af^ 
is the gauge field-strength, e the gauge coupling constant and V 
is the potential given in Eq. (H). 

This Lagrangian is invariant under the group of local U{1) trans- 
formations, 

^ e'^^'^Uix) , A^{x) ^ A^{x) + ^d,a{x) . 

The minima of the potential, (f) = r/e*^, are not invariant, the 
symmetry is spontaneously broken. Expanding as before around 
the vacuum expectation value (O|0|O) = r/, we find 

C= [d^ip - {ieA^ + id^i9){r] + ip)Wip + {leA^ + ^^^'l9){r] + if)] 
-\\^' -\f,^F^\ (11) 

where, as in the global case, rn? = Xrf. Here {} is no longer a 
physical degree of freedom. It can be absorbed by a gauge trans- 
formation. After the gauge transformation A^ A^ — {l/e)d^d 
the Lagrangian given in Eq. (|TT]) becomes 

C = {d.iff - + Im'A.A^' - \f,,F^^ + C^^t , (12) 

with m = \/~Xrj and M = \/2er]. The gauge boson "absorbs" 
the massless Goldstone boson and becomes massive. It has now 
three independent polarizations (degrees of freedom) instead of 
the original two. The phenomenon described here is called the 
'Higgs mechanism'. It works in the same way also for more com- 
plicated non-Abelian gauge theories (Yang Mills theories). 

On the classical level, what we have done here is just rewriting 
the Lagrangian density in terms of different variables. However, 
on a quantum level, particles are excitations of a vacuum state, 
a state of lowest energy, and these are clearly not described by 
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the original field but by the fields Lp and d in the global case 
and by Lp and Ai in the local case. 



The two models presented here have very close analogies in con- 
densed matter physics: 

a) The non-relativistic version of Eq. (||) is used to describe super 
fiuids where (f) is the Bose condensate (the best known example 
being super fiuid He^). 

b) The Abelian Higgs model, Eq. ([Ty) is the Landau Ginzburg 
model of super-conductivity, where represents the Cooper pair 
wave function. 

A very physical and thorough account of the problem of sponta- 
neous symmetry breaking can be found in Weinberg [|i5Y| . 

It is possible that also the scalar fields in particle physics (e.g. 
the Higgs of the standard model which is supposed to provide 
the masses of the and Z^) are not fundamental but "conden- 
sates" as in condensed matter physics. Maybe the fact that no 
fundamental scalar particle has been discovered so far has some 
deeper significance. 



2.2 Symmetry restoration at high temperature 



In particle physics like in condensed matter systems, symmetries 
which are spontaneously broken can be restored at high tem- 
peratures. The basic reason for this is that a system at finite 
temperature is not in the vacuum state which minimizes energy, 
but in a thermal state which tends to maximize entropy. We thus 
have to expand excitations of the system about a different state. 
More precisely, it is not the potential energy, but the free energy 

F = E-TS (13) 



which has to be minimized. The equilibrium value of (j) at tem- 
perature T, (0)t, is in general temperature dependent At 
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low temperature, the entropy term is unimportant. But as the 
temperature raises, low entropy becomes more and more costly 
and the state tends to raise its entropy. The field (j) becomes less 
and less ordered and thus its expectation value becomes 
smaller. Above a certain critical temperature, T > T^.^ the ex- 
pectation value vanishes. If the coupling constants are not 
extremely small, the critical temperature is of order T^. ~ 77. 

To calculate the free energy of quantum fields at finite temper- 
ature, one has to develop a perturbation theory similar to the 
T = Feynman diagrams, where ordinary Greens functions are 
replaced by thermal Greens functions. The inverse temperature, 
(3 = l/{kT) plays the role of an imaginary time component. It 
would lead us too far from the main topic of this review to intro- 
duce thermal perturbation theory, and there are excellent reviews 
on the subject available, see, e.g. [prr] , |i5Y| , ^ , |i5fc)| ,F 



Here we give a much simplified derivation of the lowest order (tree 
level) thermal correction to the effective potential . In lowest 



order the particles are non-interacting and their contributions to 
the free energy can be summed (each degree of freedom describes 
one particle), 

yeff(0,r) = F(0) + Ei^n(0,r) . (i4) 



Here V{(j)) is the zero temperature effective potential and is 
the free energy of each degree of freedom, 

= ± / ^ ln(l T exp(-e(A:)/r)) , (15) 

as known from statistical mechanics. The upper sign is valid for 
bosons and the lower one for fermions. e{k) = \Jk'^ + m^. 

For T <C nin the free energy is exponentially small. But it can 
become considerable at high temperature, T ^ m„, where we 
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obtain 



. -£r4 + ^ + C(m4) bosons 
-|gr4 + ^ + C>(m4) fermions. 



If symmetry restoration occurs at a temperature well above all 
the mass thresholds, we can approximate Ves by 

yefF(0, T) = Feff (0, T = 0) + - ^J\fT' , (17) 

Af=NB + lNF, M'= rnl + ^ E rnl . (18) 

^ bosons fermions 



Here, rrii is the formal mass given by mf = ^g^y . If the potential 
contains a 0^-term, the mass includes a term oc 0^, which leads 
to a positive quadratic term, oc T^^^. If the temperature is suffi- 
ciently high, this term overcomes the negative quadratic term in 
the Mexican hat potential and = becomes a global minimum 
of the potential. The temperature at which this happens is called 
the critical temperature. 

In the Abelian Higgs model, the critical temperature Tc, be- 
comes 

T' = (19) 



For non- Abelian 0{N) models one finds analogously \\Lb2 



l)e^ 



(20) 



The critical temperature for global symmetry breaking, i.e. with- 
out gauge field, is obtained in the limit e ^ 0. As expected, for 
~ A one finds 

Te ~ 7/ . (21) 
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Like in condensed matter systems, a phase transition is second 
order if = is a local maximum and first order if = is a 
local minimum. In the example of the Abelian Higgs model, the 
order depends on the parameters e and A of the model. 

In 0{N) models, or any other model where the vacuum manifold 
{i.e. the space of minima of the effective potential) of the bro- 
ken symmetry phase is non-trivial, minimization of the effective 
potential fixes the absolute value of but the direction, 0/|0|, 
is arbitrary. The field can vary in the vacuum manifold, given 
by the sphere S^^^ = 0{N)/0{N - 1) for 0{N) models. At low 
temperature, the free energy is minimized if the phase is constant 
(no gradient energy) but after the phase transition 4>/\4>\ will vary 
in space. The size of the patches with roughly constant direction 
is given by the correlation length ^ which is a function of time. In 
the early universe ^ is bounded by the size of the causal horizon, 

^(t) < dnit) ~ t for power law expansion. (22) 



Formally ^ diverges at the phase transition, but also our pertur- 
bative treatment is no longer valid in the vicinity of the phase 
transition since fluctuations become big. A thorough treatment 
of the physics at the phase transition is the subject of modern 
theory of critical phenomena and goes beyond the scope of this 
review. Very often, the relevant correlation length is the corre- 
lation length at the Ginsburg temperature, Tc, the temperature 
at which thermal fluctuations are comparable to the mass term. 
However, in the cosmological context there is also another scale, 
the expansion time. As the system approaches the phase tran- 
sition, it takes longer and longer to reach thermal equilibrium, 
and at some temperature, expansion is faster than the speed at 
which the system equilibrates and it falls out of thermal equilib- 
rium. It has been argued [|i(j4|] that it is the correlation length 
at this moment, somewhat before the phase transition, which is 
relevant. 
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If the phase transition is second order, the order parameter (j) 
changes continuously with time. In a first order transition, the 
state = is meta-stable (false vacuum) and the phase transi- 
tion takes place spontaneously at different positions in space and 
different temperatures T < Tc via bubble nucleation (super cool- 
ing). Thermal fluctuations and/or tunneling take the field over 
the potential barrier to the true vacuum. The bubbles of true vac- 
uum grow and eventually coalesce thereby completing the phase 
transition. 

It is interesting to note that the order of the phase transition 
is not very important in the context of defects and structure 
formation. Even though the number of defects per horizon volume 
formed at the transition does depend on the order and, especially 
on the relevant correlation length [|ib4j| , this can be compensated 
by a slight change of the phase transition temperature to obtain 
the required density of defects. 

As we have seen, a non-trivial vacuum manifold, A4 ^ {0}, in 
general implies that shortly after a phase transition the order pa- 
rameter has different values at different positions in space. Such 
non-trivial configurations are generically unstable and will even- 
tually relax to the configuration =constant, which has the low- 
est energy. Naturally, we would expect this process to happen 
with the speed of light. However, it can be slowed significantly 
for topological reasons and intermediate long lived configurations 
with well confined energy may form, these are topological defects. 
Such defects can have important consequences in cosmology. 

Several exact solutions of topological defects can be found in 
the literature, see e.g. l li-b'Zj . In the case of global defects, i.e. 
defects due to global symmetry breaking, the energy density of 
the defect is mainly due to gradient energy in the scalar field and 
is therefore not well localized in space. The scalar field gradient 
of local defects (defects due to the braking of a local, or gauge 
symmetry) is compensated by the gauge field and the energy is 
well confined to the location of the defect. To exemplify this. 
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we present the solutions for a global and a local straight cosmic 
string. 

2.3 Exact solutions for strings 
2.3.1 Global strings 

We consider a complex scalar field, (j)[x) G C, with Lagrangian 
C = \d,m-V{ct>) , V = \{\ct>\'-rff (23) 

at low temperature. The vacuum manifold is a circle of radius 
1], M. = {(t){x) G C||0| = 77}. At high temperature, T ^ rj, 
the effective potential has a single minimum at = 0. As the 
temperature drops below the critical value Tc ~ 77, a phase tran- 
sition occurs and (/) assumes a finite vacuum expectation value 
(O|0|O) 7^ which is uncorrelated at sufHciently distant points in 
physical space. If we now consider a closed curve in space 

r : [0, 1] ^ : s ^ x{s) ; a;(0) = x{l) 



it may happen that 0(a3(s)) winds around in the circle A4 ~ S^. 
We then have (f){x(s)) = r]exp(ia{s)) with a{l) = a{0) + n27r 
with n 7^ 0. Since the integer n (the winding number of the 
map r ^ : s 1-^ (J){x{s))) cannot change if we shrink the 
curve r([0, 1]) continuously, the function a{s) must be ill defined 
somewhere in the interior of F, i.e. (J) must assume the value = 
and thus have higher potential energy somewhere in the interior 

ofr([o,i]). 

If we continue this argument into the third dimension, a string of 
higher potential energy must form. The size of the region within 
which (f) leaves the vacuum manifold, the diameter of the string, 
is of the order r]~^. For topological reasons, the string cannot end. 
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It is either infinite or closed-f*! 

We now look for an exact solution of a static, infinite straight 
string along the z-direction. We make the ansatz 

(pix) = iifsipT]) exp{m(p) , (24) 



with p = and tamp = y/x^ (p is the usual polar an- 

gle. The field equation of motion then reduces to an ordinary 
differential equation for /g, 

f: + -fi - -,fs - ^fs(f! - 1) = , (25) 



where v = prj and ' = ^. A solution of this differential equation 
which satisfies the boundary conditions /s(0) = and fs{v) -^v^oo 
1 can be found numerically. It is a function of \/~Xprj. and behaves 
like 

/, - 1 - 0{l/v'^) for VXpj] > 1 ~ for VXpr] <C 1. 



The energy momentum tensor of the string is given by 



To" = T; = -^ir - \(f - If + ^,f] (26) 
T^ = Q for all other components. 

The energy per unit length of a cross-section of string out to 
radius R is 

R 

p{R) = 27r 1 Tljpdp - TTT]^ ln{VXr]R) . (27) 



The log divergence for large R results from the angular depen- 
dence of 0, the gradient energy, the last term in Eq. which 

* The only exception may occur if other defects are present. Then a string can end 
on a monopole. 
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decays only like 1/ p^. In realistic configurations an upper cutoff 
is provided by the curvature radius of the string or by its distance 
to the next string. 

Also for a single, spherically symmetric global monopoles solution 
the total energy divergies (linearly). A non-trivial results shows, 
however, that the energy needed to deform the monopole into a 
topologically trivial configuration is finite [jT[ . 



2.3.2 Local strings 

We also describe a string solution of the Abehan Higgs model, 
the Nielson-Oleson or Abrikosov vortex f[TT2 . 



The Lagrangian density is the one of scalar electrodynamics, 
Eq. (0), 



We are looking for a cylindrically symmetric, static solution of 
the field equations. For ^ oo we want the solution to approach 
a vacuum state, i.e. (f) ^ r]exp{incp) and A^^ {n/e)d^(p so that 
D^4) (the gauge field 'screens' the gradient energy). 

We insert the following ansatz into the field equations 

(f^ = lfAip)exp{in(p) (29) 

A, = ^{y/p^)aip) , Ay = -{x/p')a{p) , ^ = . (30) 
e e 

which leads to two coupled ordinary differential equations for 
and a 



d'fA ^ IdfA 
dp^ p dp p^ 



I Aid - 1; 



-jaU. 



2 

,2 2 r2, 



A 



(Pa Ida ^ ., , 
dp'^ p dp 








(31) 
(32) 
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Solutions which describe a string along the z-axis satisfy the 
asymptotics above which require 

fA{0) = a{0) = 0, and /^(p) , a{p) ^^^^ 1 . (33) 



The solution to this system of two coupled ordinary differential 
equations is easily obtained numerically. 

Asymptotically, for p — ^ oo, the a-equation reduces to the differ- 
ential equation for a modified Bessel function and we have 

a{p) c^l- pr]Ki{V2epr]), |a-l| - O (^e"^^^^) . (34) 

For large values of (3 = A/2e^, the falloff of is controlled by the 
gauge field coupling, oc (a—l)^. For (3^4 the gauge field coupling 
can be neglected at large radii p and we obtain for p — ^ oo 

fAip) - 1 - Koi^/Xpr]) - 1 - (9 (exp(-v^p77)) . (35) 

This field configuration leads to Foi = F^i = 0, hence Ei = Bi = 
B2 = ; and 

TL 

B3 = e3.,F,, = diA2 - d2Ai = —a' . (36) 



The energy per unit length of the string is 



p = 27^1^ j sds 



A 



(37) 



with s = rjp. All the terms in the integral are regular and decay 
exponentially for large s. The energy per unit length of gauge 
string is finite. The gradient energy which leads to the divergence 
for the global string is 'screened' by the gauge field. The integral is 
of the order p ~ 27r77^. This value is exact in the case (3 = \n\ = 1, 
where it can be computed analytically. In the general case with 
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|n| = 1, we have ji = 27rr]^g{/3) , where g{/3) is a slowly varying 
function of order unity. 

The thickness of a Nielson Oleson string is about 77^ and on 
length scales much larger than 77"^ we can approximate its energy 
momentum tensor by 

(T;)=;L.^(x)%)diag(l, 0,0,1) . (38) 
2.4 General remarks on topological defects 

In three spatial dimensions four different types of defects can 
form. The question whether and what kind of topological defects 
form during a symmetry breaking phase transition is determined 
by the topology of the vacuum manifold A4: 

• If Ai is disconnected, domain walls from. Example: if the 
symmetry (J) ^ —(J) for a real scalar field is spontaneously bro- 
ken, A4 = {—77,77}. Domain walls form when discrete symme- 
tries are broken. Discrete symmetries are not continuous and 
therefore cannot be gauged. Hence domain walls are always 
global defects. 

• If there exist loops in A4 which cannot be continuously shrunk 
into a point, Ai is not simply connected, strings form. 
Example: if U{1) is completely broken by a complex scalar 
field, = S^, see previous subsection. 

• If A4 contains non-contractible spheres, monopoles form. 
Example: if 0(3) is broken to 0(2) by a three component scalar 
field, M = S^ = 0(3)/0(2). 

• If A4 contains non-contractible 3-spheres, textures form. 
Example: if 0(4) is broken to 0(3) by a four component scalar 
field, 7W = S3 = 0(4)/0(3). 

These topological properties of Ai are best described by the 
homotopy groups, 7Tn{Ai)- The group 7i^{A4) is relevant for the 
existence of textures, 7T2{M) decides about monopoles, 7Ti{M.) is 
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relevant for strings and 7Tq{A4) for domain walls []S^]. If a symme- 
try group G is spontaneously broken to a subgroup H G the 
vacuum manifold is in general equivalent to the quotient space, 
A4 ~ G/H. In the monopole example above, we have G = 0(3) 
and H = 0{2). The vacuum manifold is ~ ~ 0(3)/0(2). 



2. 5 Defect formation and evolution in the expanding universe 



Our universe which is to a good approximation an expanding 
Friedman universe was much denser and hotter in the past. Dur- 
ing expansion the universe may cool through a certain critical 
temperature T^. at which a symmetry G is spontaneously broken 
down to H C G. If A4 = G/H is topologically non-trivial, topo- 
logical defects can form during the phase transition. This scenario 
is called the Kibble mechanism . We apply the Kibble mecha- 
nism to estimate the energy density in defects from phase transi- 
tions with different vacuum manifolds at a given temperature Tc. 
Consider a field theory with symmetry group G and Higgs-field 
(f) with a self- interaction potential V{(j)). For illustration we use 
G C, G = U{1) and 

= iA(# - tjY- (39) 

At finite temperature, the free energy is of the form 

^r(0) = ^r200 + H(0), (40) 



where ^4 is a real constant given by combinations of A and other 
coupling constants {e.g. gauge couplings, Yukawa couplings). 
The sign of A depends on the number of fermions. We assume 
^4 > 0, i.e. that there are only few fermions and sufficiently 
small Yukawa- couplings. Then, from Eqs. ( 13^1 ) and (|5y), we see 
that the effective masses of the field at temperature T and Tc 
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are 



m\T) = "\4'(0 = 0) = AT'^ - Xti\ 
m(Te) = 0, Te = r/ (^-j , for A - ^ - 1, - r]. 



At r = Tc, this field theory undergoes a second order phase 
transition: the equihbrium point = becomes unstable for 
T < Tc {rn? becomes negative) and assumes a non-vanishing 
vacuum expectation value. 

For another form of Vt-, the equilibrium = at T = can be 
meta-stable so that the phase transition is of first order. Hence, 
to decide whether the transition is of first or second order, it is 
important that we can rely on the form of the effective potential 
Vr(0) which is obtained by perturbation theory or by numerical 
lattice calculations. This is in general a difficult problem. For the 
electro- weak theory, e.g. , it was discovered only recently that 
the electroweak transition is probably not a real phase transition 
but just a continuous cross-over . 



The correlations of the field are described by the thermal 
Greens functions: 

G{\x-x'\) = {ct)\x,mx',t)). 



For massive particles at T > Tc where cjf = /c^ + m^{T) one can 
write 



1 e*(^'^') 



G(|.-.'|) = 2/^^^ + G», 
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where are the zero temperature contributions. For T ^ Tc 
such that m{T) <C T we have, with r = |a3 — cc'l: 



G{r) ~ 



rV6, forr<Cl/r 
exp[-m(r)r]r/(27rr), for r > l/T 



For T Tc, m{T) and therefore, at Tc G ~ for large 

r. The correlation length ^ for the phase transition (of 2"^ order) 

is defined as ^ := t^t. This definition is different from the defini- 
te 

tion used in solid state physics. There one defines the correlation 
length as the length above which the correlation decreases expo- 
nentially. In this sense, the correlation length would be infinite 
at Tc {vc = l/m(r)). In cosmology the correlation length cannot 
diverge because of causality. It is bounded from above by the dis- 
tance a photon can travel during the age of the universe until tc- 
This distance is (for non-inflationary expansion) given by 

Ih = a{r)T = tc: tc ■= / a{r) dr. 



Hence another meaningful definition of the correlation length 
would he ^ := Ih = tc. Often also the correlation length at the 
Ginzburg temperature or at the temperature (before the phase 
transition) at which the system drops out of thermal equilib- 
rium |i(34|] is chosen. For the following it is not important which 



of the above definitions we use. We only require that ^ < Ih- 
We now suppose that directly after the phase transition, the vac- 
uum expectation value (0) takes arbitrary uncorrelated values in 
points with distance r > ^, but stays continuous (finite gradient 
energy!). If 7Tn{A4) is non-trivial for n < 3, the map 

(0) : S"^>(, (0(cc)) 



for a large enough n-sphere S" in physical space, may represent a 
non-trivial element of 7r„(A1). Then 0(S") cannot be contracted 
continuously to a point on A4 and, somewhere inside S*^, (0) has 



28 



to leave the vacuum manifold, {4>){p) ^ A4. These positions of 
higher potential energy are topological defects. The type of defect 
formed depends on the order n of the non-trivial homotopy group: 

• n = 0: 2-dimensional defects, domain walls, d = 3 

• n = 1: line-like defects, cosmic strings, d = 2 

• n = 2: point-like defects, monopoles, d = 1 

• n = 3: event-like defects, texture, d = 0. 

Here d is the spacetime dimension of the defect, d = 4 — 1 — n. 
If a vacuum manifold has several nontrivial homotopy groups 
with n < 3, generically only the lowest n defects survive and the 
higher order defects are instable. As an example, in the isotropic 
to nematic phase transition of liquid crystals 0{3) is broken 
to 0(2) X Z2 leading to M = 0(3)/(0(2) x Z2 = 87^2- This 
allows for texture, monopoles and strings, but textures decay 
into monopole anti-monopole pairs and monopole/anti-monopole 
pairs are connected by strings and attract each other until they 
annihilate. Only strings scale [j2^ . 



The case of texture, n = 3 can be described in this context only if 
either the universe is closed and physical space is a three sphere 

of if (j) is asymptotically parallel, i.e. (I){x^t) '"^^ — >^ 0o- Then the 
points |cc| ^ 00 can be identified in all directions and we can re- 
gard as a map from U {00} = to and ask whether this 
map is topologically trivial or not. In the cosmological context 
this concept violates causahty. However, the texture case allows 
for a texture winding number density whose integral over all of 
space only takes integer values if is asymptotically constant (or 
space is a three sphere). The integral of the winding number den- 
sity over a region of space tells us whether the field configuration 
inside contains textures. 

According to the above description of the process of defect for- 
mation after a cosmological phase transitions, called the Kibble 
mechanism [^j, we typically expect on the order of one defect 



per horizon volume. Simulations and analytical arguments show 
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that the actual number is somewhat larger for cosmic strings and 
significantly smaller for texture. 



If defects are local, the scalar field gradiants are compensated 
by the gauge field and they do not interact at large distance 
other than gravitationally in the simplest model, where no mass- 
less charged particles 'five' inside the defect. An exeption to 
this are superconducting cosmic string [ |ib9|j . For example local 



monopoles do not annihilate once they are formed and their den- 
sity just scales with the expansion of the Universe, like 
Since they are non-relativistic, m > T, their energy density scales 
the same way and they soon dominate the total energy density 
of the universe. Every simple GUT group produces monopoles 
when it breaks down to the standard model symmetry group, 
SU (3) X SU (2) X [/(I). The observed absence of monopoles there- 
fore represents a serious problem for the unification of standard 
cosmology with grand unified theories Local texture, on the 
contrary, soon thin out and do not induce sufficiently strong 
perturbations to generate structure in the universe. Only local 
strings scale, i.e. contribute a constant small fraction to the en- 
ergy density of the universe, and are therefore possible candi- 
dates of topological defects for structure formation. If the group 
7ri(A^) is non-Abelian (tti is the only homotopy group which can 
be non-Abelian), the cosmic string network becomes 'frustrated' 
and does not scale. Such a low energy frustrated string network 
has been proposed as candidate for the cosmic dark energy 



The situation is different for global defects. There, the main con- 
tribution to the energy density comes from the Higgs field and 
scales as 1/^^, like the background energy density in the universe 
(up to logarithmic corrections for global strings). The only excep- 
tion are domain walls which are forbidden, since they soon come 
to dominate, leading to a very inhomogeneous universe. Recently, 
however 'soft domain walls' 0], i.e. domain walls forming at a late 
time phase transition, have been studied. 
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3 Theoretical Framework 



3.1 Linear cosmological perturbations with seeds 



A basic tool for cosmic structure formation is linear cosmological 
perturbation theory. The fact that CMB anisotropics are small 
shows that at least initially also perturbations in the matter den- 
sity have been much smaller than unity and therefore they may 
be treated within linear perturbation theory. It is generally as- 
sumed (an assumption which is supported by several observa- 
tional facts, see e.g. [ |i5Yi| ) that perturbations are still linear on 



scales above about 10/i^^Mpc. On smaller scales non-linear N- 
body simulations are needed to compute the evolution of density 
fluctuations. 

The principal difference in perturbation theory in models with 
topological defects as compared to the more familiar inflation- 
ary models, is the fact that here cosmic perturbation equations 
are not homogeneous. The perturbations are induced by 'seeds' 
which are not present in the background energy momentum ten- 
sor.. The defect energy momentum tensor enters in the pertur- 
bation equation as 'source' or 'seed' term on the right hand side, 
but the defects themselves evolve according to the background 
space-time. Perturbations in the defect evolution are of second 
order. (This procedure has sometimes also been termed the 'stiff 
approximation' [|i5U|] , but it is actually nothing else than consis- 



tent linear perturbation theory.). 

Gauge- invariant perturbation equations for cosmological models 
with seeds have been derived in Refs. [^,^. Here we follow the 



notation and use the results presented in Ref. [j52j . Definitions of 
all the gauge-invariant perturbation variables used in this Review 
in terms of perturbations of the metric, the energy momentum 
tensor and the photon and neutrino brightness are given in Ap- 
pendix A for completeness. 
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We consider a background universe with density parameter Qq = 
Q„i + r^A = 1, consisting of photons, cold dark matter (CDM), 
baryons and neutrinos. At very early times z ^ z^ec ~ 1100, 
photons and baryons form a perfectly coupled ideal fluid. As time 
evolves, and as the electron density drops due to recombination 
of primordial helium and hydrogen, Compton scattering becomes 
less frequent and higher moments in the photon distribution de- 
velop. This process has to be described by a Boltzmann equation. 
Long after recombination, free electrons are so sparse that the col- 
lision term can be neglected, and photons evolve according to the 
coUisionless Boltzmann or Liouville equation. During the epoch 
of interest here, neutrinos are always coUisionless and thus obey 
the Liouville equation. 

In the next section, we parameterize in a completely general way 
the degrees of freedom of the seed energy momentum tensor. 
Section [3T3| is devoted to the perturbation of Einstein's equations 
and the fluid equations of motion. Next we treat the evolution 
of CMB photons by the Boltzmann perturbation equation, in- 
cluding polarization. The detailed derivations as well as the ex- 
pressions for the CMB anisotropy and polarization power spectra 
are given in Appendix p. We finally explain how to compute the 
power spectra of density fluctuations, CMB anisotropics and pe- 
culiar velocities by means of the derived perturbation equations 
and the unequal time correlators of the seed energy momentum 
tensor which are obtained by numerical simulations. 



3.2 The seed energy momentum tensor 



Since the energy momentum tensor of the seeds, G^j,, does not 
contribute to the background Friedman universe, it is gauge in- 
variant by itself according to the Stewart- Walker Lemma [|i4U| . 



0^1, can be calculated by solving the matter equations for the 
seeds in the Friedman background geometry. Since G^jy has no 
background component it satisfies the unperturbed "conserva- 
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tion" equations. We decompose Q^^, into scalar, vector and ten- 
sor contributions. They decouple within linear perturbation the- 
ory and it is thus possible to write the equations for each of 
these contributions separately. As always (unless noted other- 
wise), we work in Fourier space, k is the comoving wave number 
and k = \k\. We parameterize the scalar (S) vector {V) and 
tensor (T) contributions to G^j^ in the form 



^jo - 



1 



(/p + f^)6ji - kjkiU 



= iM^l(kjwi''^ + k 



(vr) 



(41) 
(42) 

(43) 

(44) 

(45) 

(46) 



Here M denotes a typical mass scale of the seeds. In the case 
of topological defects we set M = 77, where 7/ is the symmetry 
breaking scale The vectors w'^'"^ and w^'^^ are transverse and 



r/j^'' is a transverse traceless tensor. 



k ■ w^""^ = k ■ w^""^ = kW^J^ = ^ ' = 



From the full energy momentum tensor G^j, which contains scalar, 
vector and tensor contributions, the scalar parts fy and fj^ of a 
Fourier mode are given by 

ik^eoj = -k^M^fy , (47) 
- k'k\e,, - ^(5,,(5^^GfcO = ^k'M^U . (48) 

On the other hand fy and fj^ are also determined in terms of 
fp = Goo/M^ and fp = G^j/(3M^) by energy and momentum 
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conservation, 



fp + k'f, + -{fp + 3fp) = , (49) 
a 

A + -f, + h'U = . (50) 



Once fy is known it is easy to extract 



M'^wf^ = Soj - ikjM^fy . (51) 



boT wl we use 



^k^ie, - e<f ) = . (52) 



Again, wj^^ can also be obtained in terms of W;""* by means of 
momentum conservation, 

+ + -k^w^'''^ = . (53) 

a 2 



Finally, 



k' 



{fp+Jfn)S^J-k,kJU 

l:{kiwf^ + kjW^^^) 



(54) 



The geometry perturbations induced by the seeds are charac- 
terized by the Bardeen potentials, $s and ^g, for scalar per- 
turbations, by the potential for the shear of the extrinsic cur- 
vature, for vector perturbations and by the gravitational 

is) 

wave amplitude, H\j\ for tensor perturbations. Detailed defini- 
tions of these variables and their geometrical interpretation are 
given in Ref. [ji^] (see also Appendix A). Einstein's equations 
for an unperturbed cosmic background fluid with seeds relate 
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the seed perturbations of the geometry to the energy momentum 
tensor of the seeds. Defining the dimensionless small parameter 



e = AtiGM'^ , (55) 
we obtain to first order in e 



k'^s = e{fp + ?>-h) (56) 
d 

$, + *, = -2eA (57) 

-A;2Ef) = 4e4"^ (58) 

H\f + 2-H\f + k'Hlf = 2eT^f . (59) 



Eqs. (^B]) to ( |5D| ) would determine the geometric perturbations 



if the cosmic fluid were perfectly unperturbed. In a realistic sit- 
uation, however, we have to add the fluid perturbations which 
are defined in the next subsection. Only the total geometrical 
perturbations are determined via Einstein's equations. In this 
sense, Eqs. (|5^) to (^ should be regarded as definitions for 
,E(^) and H\f. 

A description of the numerical calculation of the energy momen- 
tum tensor of the seeds for global defects and cosmic strings is 
given in Chapter 



3. 3 Einstein's equations and the fluid equations 



3. 3. 1 Scalar perturbations 

Scalar perturbations of the geometry have two degrees of freedom 
which can be cast in terms of the gauge-invariant Bardeen poten- 
tials, ^ and $ [1§,|9T|. For Newtonian forms of matter, ^ = — $ 



is nothing else than the Newtonian gravitational potential. For 
matter with significant anisotropic stresses, ^ and — $ differ. In 
geometrical terms, the former represents the lapse function of the 
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zero-shear hyper-surfaces while the latter is a measure of their 3- 
curvature [^. In the presence of seeds, the Bardeen potentials 
are given by 



^ = *. + ^m, (60) 
$ = $. + $m, (61) 

where the indices s,m refer to contributions from a source (the 
seed) and the cosmic fluid respectively. The seed Bardeen poten- 
tials are given in Eqs. (^ and (0). 

To describe the scalar perturbations of the energy momentum 
tensor of a given matter component, we use the gauge invari- 
ant variables Dg for density fluctuations, corresponding to the 
usual density fluctuation in the 'flat gauge', V, for the potential 
of peculiar velocity fluctuations, corresponding to the usual ve- 
locity potential in the longitudinal gauge and 11, a potential for 
anisotropic stresses (which vanishes for CDM and baryons). A 
definition of these variables in terms of the components of the 
energy momentum tensor of the fluids and the metric perturba- 
tions can be found in Refs. [^Tj or and in Appendix A. 



Subscripts and superscripts ^, c, 6 or denote the photons, CDM, 
baryons and neutrinos respectively. 

Einstein's equations yield the following relation for the matter 
part of the Bardeen potentials 



+3pc+3pb +4/?^$ +3-k-^{^p^V^ +pcVc +pbVb +|/9^K}I62) 

*m = -«l>m ^ {pjU^ + p,,U^) . (63) 

Note the appearance of $ = $s + on the r.h.s. of Eq. (|62|). Us- 
ing the decompositions (^,0) we can solve for $ and ^ in terms 
of the fluid variables and the seeds. With the help of Friedman's 
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equation, Eqs. ( |62|) and ( |63|) can then be written in the form 



$ = o [x^D'^^^ + Xe^l") + XD^^) 

3 la j /C +4x^+3a;c+3a;b+4x^ 



^ a 
2 



a 

-2 



+4x,K) + 3^:^ (^) ^s] (64) 
^ = - 2e/^ - Q A;"2(x^n^ + x^H^) . (65) 



Here we have normahzed the scale factor such that a = 1 today. 
The density parameters Q, always represent the values of the 
corresponding density parameter today (Here , stands for c , 7 , 6 
or J,.). To avoid any confusion, we have introduced the variables 
X, for the time dependent density parameters. 



n 



X^^y 



Xc.b 



fl^ + flea + + r^i/ + ^Att"^ 
^cbd 

fl^ + flcO' + ^bd + + Qa^^ 



(66) 
(67) 



The fluid variables for photons and neutrinos are obtained by 
integration over directions of the scalar brightness perturbations, 
which we denote by A4s{T,k,n) and Ms{T,k,ri) respectively. 
They are given in Appendix 0. 

The evolution of CDM perturbations is determined by energy 
and momentum conservation, 

D'g''^ + kVc = , Vc+ (^) Vc = k^ . (68) 



During the very tight coupling regime, z ^ Zdec-, we may neglect 
the baryon contribution in the energy momentum conservation 
of the baryon-photon plasma. We then have 
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+ ^kV, = , V,- k^Dp^ = k{^ - , (69) 
and Df^ = ^D^;\ V, = V,. (70) 

The conservation equations for neutrinos are not very useful, 
since they involve anisotropic stresses and thus do not close. At 
the temperatures of interest to us, T <C IMeV, neutrinos have to 
be evolved by means of the Liouville equation which we discuss 
in the next section. 

Once the baryon contribution to the baryon-photon fluid be- 
comes non-negligible, and the imperfect coupling of photons and 
baryons has to be taken into account (for a 1% accuracy of the re- 
sults, the redshift corresponding to this epoch is around z ~ 10''), 
we evolve also the photons with a Boltzmann equation. The equa- 
tion of motion for the baryons is then 

Df^ + kVi, = 0, (71) 

H+(«)vi = fe*-^^5^[v,-yd. (72) 

\aj 6\li) 

The last term in Eq. ([72|) represents the photon drag force in- 
duced by non-relativistic Compton scattering, cry is the Thomson 
cross section, and denotes the number density of free electrons. 
The scale factor a enters since our derivative is taken w.r.t confor- 
mal time. At very early times, when arn^a ^ 1/r, the 'Thomson 
drag' just forces V5 = Vr^, which together with Eqs. ( 1691 ) and ( [711 ) 
implies the first eqn. of ([701). 

An interesting phenomenon often called 'compensation' can be 
important on super horizon scales, /cr <C 1. If we neglect aniso- 
tropic stresses of photons and neutrinos and take into account 
that 0{V) = C(A;r*) and 0{Dg) = O(^) for kr < 1, Eqs. (0) 
and (^) lead to 

0($) = O {{krf^, - 2eA) . (73) 
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Hence, if anisotropic stresses are relatively small, e/^r <C $s, the 
resulting gravitational potential on super horizon scales is much 
smaller than the one induced by the seeds alone. One must be 
very careful not to over interpret this 'compensation' which is 
not strictly related to causality, but is due to the initial condi- 
tion Dg , V -^T^o 0. A thorough discussion of this issue is found 
in Refs. [^ , |23| , p^] . As we shall see, for textures $s and efj^ are 



actually of the same order. Therefore Eq. ( |73|) does not lead to 
compensation, but it indicates that CMB anisotropics on very 
large scales (Sachs- Wolfe effect) are dominated by the amplitude 
of seed anisotropic stresses. Nevertheless, for purely scalar or co- 
herent perturbations, as we shall see /^r oc (/cr)^$s and hence 
'compensation' is important. 

The quantities which we want to calculate and compare with ob- 
servations are the CDM density power spectrum and the peculiar 
velocity power spectrum today 

P(/c) = (|D(^)(A:,ro)|2) and P„(A:) = ro)P) . (74) 

Here (• • •) denotes an ensemble average over models. Note that 
even though Dg and V are gauge invariant quantities which do 
not agree with, e.g., the corresponding quantities in synchronous 
gauge, this difference is very small on sub-horizon scales (of order 
1 / kr) and can thus be ignored. 

On sub-horizon scales the seeds decay, and CDM perturbations 
evolve freely. We then have, like in inflationary models [|ii5|] , 

P,{k) = Hl^l'P{k)k~' . (75) 



3. 3. 2 Vector perturbations 

Vector perturbations of the geometry have two degrees of free- 
dom which can be cast in a divergence free vector field. A gauge- 
invariant quantity describing vector perturbations of the geome- 
try is E, a vector potential for the shear tensor of the {r =const.} 
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hypersurfaces. Like for scalar perturbations, we split S into a 
source term coming from the seeds given in the previous section, 
and a part due to the vector perturbations in the fluid, 

E = + . (76) 



The perturbation of Einstein's equation for E^ is 

/C^E^ = 6 Q i^X^LD^ + XcUJc + XbUJb + ^X„UJ„] . (77) 



Here uj, is the fluid vorticity which generates the vector type 
shear of the equal time hyper-surfaces (see Appendix A). By def- 
inition, vector perturbations are transverse, 

T.-k = T>m-k = T>s-k = u),-k = 0. (78) 



It is interesting to note that vector perturbations in the geom- 
etry do not induce any vector perturbations in the CDM (up 
to unphysical gauge modes), since no geometric terms enter the 
momentum conservation for CDM vorticity, 

a 
a 



hence we may simply set Uc = 0. This is also the case for the 
tightly coupled baryon radiation plasma. But as soon as higher 
moments in the photon distribution build up, they feel the vector 
perturbations in the geometry (see next section) and transfer it 
onto the baryons via the photon drag force. 



ujb^[-jujb = — [uj^-ujb]. (79) 



The photon vorticity is derived via the integral over the vector 
type photon brightness perturbation, A4v, using 

^7 = v7 _ (80) 
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and 

= — [ nM^^^dn , (81) 

where the integral is over photon directions, n (see Appendix IB]). 
The vector equations of motion for photons and neutrinos are 
discussed in the next section. 

3. 3. 3 Tensor perturbations 

Metric perturbations also have two tensorial degrees of freedom, 
gravity waves, which are represented by the two helicity states 
of a transverse traceless tensor (see Appendix A). As before, we 
split the geometry perturbation into a part induced by the seeds 
and a part due to the matter fluids, 

H„ = ff<;' + Hf . (82) 



The only matter perturbations which generate gravity waves are 
tensor type anisotropic stresses which are present in the pho- 
ton and neutrino fluids. Numerically one finds that the effect 
of anisotropic stresses of photons and neutrinos contributes less 
than 1% to the final result [^, and hence may be neglected by 
setting h\J^^ = 0. 



3.4 Boltzmann equation, polarization and CMB power spectra 



When particle interactions are less frequent, the fluid approxima- 
tion is not sufficient, and we have to describe the given particle 
species by a Boltzmann equation, in order to take into account 
phenomena like collisional and directional dispersion. In the case 
of massless particles like massless neutrini or photons, the Boltz- 
mann equation can be integrated over energy, and we obtain an 
equation for the brightness perturbation which depends only on 
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momentum directions [^. As before, we split the brightness per- 
turbation into a scalar, vector and tensor component, and we 
discuss the perturbation equation of each of them separately, ^ 



M = Ms + Mv + Mt 



(83) 



The function A4 depends on the wave vector fc, the photon di- 
rection n and conformal time r. Linear polarization of photons 
induced by Compton scattering is described by the Stokes pa- 
rameter Q and depending on the same variables. An explicit 
derivation of the Boltzmann equation including polarization is 
presented in Appendix p. Here we just repeat the necessary def- 
initions and results. 

The brightness anisotropy A4 and the non-vanishing Stokes pa- 
rameters Q and U can be expanded as 

X(r,fc,n) = E E M^r\r,k)oGT{n), (84) 

e TO=-2 

g(r,/c,n)±zt/(r,fc,n) = E E (^f^ ± z^f ^)2G',-(n).(85) 

e TO=-2 

The spin weighted spherical harmonics are described in Ap- 
pendix 0. Up to a normalization constant, the qG^ coincide with 
the usual spherical harmonics. The coefficients m = 0, m = ±1 
and m = ±2 describe the scalar (S), vector {V) and tensor (T) 
components respectively. The Boltzmann equation for the coeffi- 
cients Xf^^ is given by 



k 



2i-l 

(m) q{m) 



2£ + 3 
{I > m) 



(86) 



**We could in principle add higher spin components to the distribution functions. 
But they are not seeded by gravity and since photons (and neutrinos) interact at 
high enough temperatures, they are also absent in the initial conditions. 
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(88) 



where we set 



((0) 



5f 



nedraC^o), 5}'^ = ne(Jra4a;6, (89) 

5^1^ = nedTflC^i) + 4E, Si^^ = UearaC^^'^ + 4i7 



and C^"*) = i7t[A^2'"^ — VQE2^^]. The couphng coefficients are 



i{m) 



sK,£ — 



\ 



(£2-m2)(£2-g2; 



In Appendix 0) we express the fluid variables in terms of integrals 
of the photon brightness over directions. 

The CMB temperature and polarization power spectra are given 
in terms of the expansion coefficients A4^P\ E^^^ and Bg"^^ as 

{2£ + Ifcf^^"^^ = 1^ / /c^d/cXf V/™)* , (90) 



where Um = 1 for m = and = 2 for m = 1, 2, accounting 
for the number of modes. Since B is parity odd, the only non- 
vanishing cross correlation spectrum is C^^. 
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3. 5 Neutrinos 



Analogously to the photon brightness perturbation it is useful to 
introduce the neutrino one as well, which we will call A/", 



A/- = A/-<*>+A/-"')+A/-<^>. 



(91) 



Since the neutrinos are collisionless during the entire epoch un- 
der consideration, they satisfy the collisionless equations which 
are obtained from the Boltzmann equation for the intensity by 
setting gt = 0. For simplicity, we expand M not in terms of the 
functions QG'f^\ but use the more basic approach with Legendre 
polynomials. 



J\fiv) = ^i_^2 ]^j^iv) 0+a4^^ ill) sin 

J^m = {l-n^) [a^''' cos(2(/.) + ATf * sin(2(^) 
JVS = E(-«)'(2^+l)^^r+,x)^'(M). 



(92) 

(93) 
(94) 

(95) 
(96) 



The Liouville equation for the coefficients i^jf '^'"^^ then becomes 



AS) 


k 






-((+1)4% 






k 






-ie + i)uZ{ 






k 






2£ + l 



= -/c(*-$)^,,i (97) 
= 4k6i,ii: (98) 
= 4i^e . (99) 
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3.6 Computing power spectra in seed models 



The generation of the seeds, e. 5'. topological defects during a 
symmetry breaking phase transition, is an inherently random 
process. The exact seed distribution in our universe is just one 
realization and cannot be predicted. Only statistical properties, 
expectation values, can be calculated. Yet the source functions 
to the Boltzmann equation are elements of the seed energy mo- 
mentum tensor, not their expectation values. 

In principle one could calculate the induced random variables 
Df\k,TQ), Vcik.To), M'r\k,To) etc for 100 to 1000 realizations 
of a given model and determine the expectation values P{k), 
Pv{k) and Ci by averaging. This procedure has been adapted in 
Ref. [0] for a seed energy momentum tensor modeled by a few 
random parameters and in Ref. [j^ where the CMB anisotropics 



on large scales have been determined in cc-space by direct line of 
sight integration and averaging over several observer positions. 

In a more realistic calculation of topological defects, where the 
seed energy momentum tensor comes entirely from numerical sim- 
ulations, this procedure is not feasible. The first and most impor- 
tant bottleneck is the dynamical range of the simulation which 
is about 40 in the largest (512)^ simulation which have been per- 
formed [[n7| ,|47|] . They need about 1 to 2 Gbyte of RAM and run 
in about one hour CPU time on a modern work station or PC. To 
determine the C^s for 2 < £ < 1000 we need a dynamical range of 
about 10,000 in /c-space. This means /Cmax/^min ~ lO'OOO, where 
/Cmax and /Cniin circ the maximum and minimum wave numbers 
which contribute to the C^'s to achieve and accuracy of about 
10%. A dynamical range of 10,000 requires at least a (100,000)^ 
simulations which needs about 10,000 Terabytes RAM! Corre- 
spondingly the CPU time required for such a simulation is about 
1000 years. 

With brute force, this problem is thus not tractable with present 
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or near future computing capabilities. But there are a series of 
theoretical observations which reduce the problem to a feasible 
one: 

As we have seen, for each wave vector k given, we have to solve 
a system of linear perturbation equations with random sources, 

VX = S . (100) 

Here X> is a time dependent linear differential operator, X is 
the vector of the matter perturbation variables specified in the 
previous subsections (photons, CDM, baryons and neutrini; total 
length up to 4000), and S is the random source term, consisting 
of linear combinations of the seed energy momentum tensor. 

For given initial conditions, this equation can be solved by means 
of a Green's function (kernel), ^(r, r'), in the form 

To 

X^(ro, k) = J drgj„i{ro, r, k)S^^{T, k) . (101) 



We want to compute power spectra or, more generally, quadratic 
expectation values of the form 

(X,-(ro,fc)X;(ro,fcO) , 
which, according to Eq. (|i01|) are given by 



(X,(ro,fc)X;(ro,fc')) = 

JdrGjUr. k) Jdr'giir', k') x (<S^(r, fc)<S:(r', k')) . (102) 

The only information about the source random variable which 
we really need in order to compute power spectra are therefore 
the unequal time two point correlators 

{Srn{r,k)S:{T\k')) . (103) 
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This nearly trivial fact has been introduced by Hindmarsh [[70 
and exploited by many workers in the field. For example in Ref. 



where decoherence of models with seeds has been discovered, 
and later in Refs. [|i 1 Y| ,^, P5| , ^ , ^] and others. The eigenvector 



method discussed here has been introduced in [|i4;j|] . (The CMB 
anisotropy spectrum from cosmic texture shown in this paper is, 
however, incorrect.) 

To solve the enormous problem of dynamical range, one then uses 
causality, statistical isotropy and 'scaling'. 

Seeds are called 'scaling' if their correlation functions C^^^px de- 
fined by 

e,,{k,r) = M^9,,{k,r) , (104) 
(V(fc, r)e;,{k\ r')) = Cp^pxik, r, r')6{k - k') (105) 

are scale free; i.e. the only dimensional parameters in C^j^p\ are 
the variables r, t' and k themselves. The ^-function in fc-space 
is a simple consequence of statistical homogeneity. Up to a cer- 
tain number of dimensionless functions F„ oi z = k\/rr' and 
r = t/t', the correlation functions are then determined by the 
requirement of statistical isotropy, symmetries and by their di- 
mension. Causality requires the functions F„ to be analytic in 
. A more detailed investigation of these arguments and their 
consequences is given in Chapter |^. There we also show that sta- 
tistical isotropy and energy momentum conservation reduce the 
correlators ( |105| ) to five such functions Fi to F5. 

In cosmic string simulations, energy and momentum are not con- 
served. Cosmic string loops oscillate and emit gravitational waves 
(see Refs. [|151| , CT] ). They lose their energy by radiation of grav- 



itational waves and, in 'cusps' or tiny wiggles into massive par- 
ticles f|lb2|]. In this case 14 functions of z'^ and r are needed to 



describe the unequal time correlators 



Since analytic functions generically are constant for small argu- 
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ments <^ 1, F„(0, r) actually determines F„ for all values of 
k with z = k^frT' ~ 0.5. Furthermore, the correlation func- 
tions decay inside the horizon and we can safely set them to 
zero for 2; ~ 40 where they have decayed by about two orders of 
magnitude (see Figs. || to in Chapter |]). Making use of these 
generic properties of the correlators, we have reduced the dynam- 
ical range needed for our computation to about 40, which can be 
attained with (512)^ simulations feasible on present computers. 

Clearly, all correlations between scalar and vector, scalar and 
tensor as well as vector and tensor perturbations have to vanish. 

The source correlation matrix Cy^ypuih^ r, r') can be considered as 
kernel of a positive hermitian operator in the variables x = kr = 
zr^/'^ and x' = kr' = z/r^^'^^ which can be diagonalized: 

C{x,x') = Y.Kvn{x)vl{x') (106) 

n 



(the variable k and the space time indices are supressed in this 
and the following expressions). The series [vn) is an orthonormal 
series of eigenvectors of the operator C (ordered according to the 
amplitude of the corresponding eigenvalue) for a given weight 
function w. We then havej^^^ 

j C{x,x')v^''\x')w{x')dx' = XnV^'^^x) . (107) 



The eigenvectors and eigenvalues depend on the weight function 
w which can be chosen to optimize the convergence speed of the 
sum ( llObi) . For 0{N) models, scalar perturbations typically need 



20 eigenvectors whereas vector and tensor perturbations need five 
to ten eigenvectors for an accuracy of a few percent (see Fig. § 
in Chapter |5|). 



* '^itkie the assumption that the operator C is trace-class enters. This hypothesis is 
verified numerically by the fast convergence of the sum (|106|) . 
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Inserting the expansion ( |iut)|) in Eq. (|iu:^|) , leads to 



{X,{k,To)X*{k,ro)) =Y.KXt'\kro)X)"^\krQ) , (108) 



where xj;^\To) is the solution of Eq. (llUUl ) with deterministic 



source term v. 



in) 



(109) 



For the CMB anisotropy spectrum this gives 

C, = Y. A!."'cf + £ Af'cf + £ A,f . (110) 



C * is the CMB anisotropy induced by the deterministic source 



v^*'^' , and n, is the number of eigenvalues which have to be con- 
sidered to achieve good accuracy. Here we have also used that 
the unequal time correlation matrix contains uncorrelated scalar, 
vector and tensor blocks. 



Instead of averaging over random solutions of Eq. (|iUi| ), we can 
thus integrate Eq. ( |iUl| ) with the deterministic source term v^'^^ 
and sum up the resulting power spectra. The computational re- 
quirement for the determination of the power spectra of one seed 
model with given source term is thus on the order of ns + nv + UT 
inflationary models. This eigenvector method has first been ap- 
plied in Ref . [|ll Y| . 



This completes the formal developments needed to compute struc- 
ture formation with defects. In the next chapter we discuss the 
numerical simulations which have been performed to obtain the 
unequal time correlators of the seed energy momentum tensor. 
These are then diagonalized and the eigenfunctions are entered 
as sources in the system of linear equations derived in this chap- 
ter. The results from this procedure are described in Chapter 
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4 Numerical Implementation 



In the previous chapter we have learned that the only input 
needed for the computation of power spectra and other two point 
correlation functions are the unequal time correlators of the de- 
fect energy momentum tensor. In this chapter we discuss how 
they are obtained in practice. Their general structure will be 
analysed in Chapter 

4-1 Global defects 

4.1.1 The a model approximation 

We consider a spontaneously broken scalar field with O(A^) sym- 
metry. If we are not interested in the microscopical structure of 
the field in the vicinity of the core but only in its behaviour on 
large scales, we can force the field to stay on the vacuum manifold 
with a Lagrange multiplier A and drop the potential. The bulk 
part of the energy of global strings, global monopoles and global 
texture is contained in the field gradient at large distances of the 
core and is not affected by this approximation. 

c = d^4)-d^4) + \{(t? -rf) . (Ill) 

Varying the action with respect to and A leads to 

□0 + A0 = O, 02-r/2 = O. (112) 

Multiplying the first equation with and using cjp' = rf yields 
A = (0- □0)/r/2 hence 00 - (0 • □0)0/7/^ = . (113) 

Applying d^d^^ on 0^ — rf we find (0 • D^) = —{d^cf) ■ d^cf)) so that 

□0+(a^(/)-6"^0)0/7/2 = O . 
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Setting = 0/?7, we finally obtain the equation of motion 

nf]+{df,p-d^f3)f3 = (114) 

for the field (3 G S*^, where n = 1,2 and 3 for global strings, 
monopoles and texture respectively. Eq. ( |114| ) is the equation of 



motion for the non-linear cr-model for a scalar field on S 



n 



4-. 1.2 The energy momentum tensor of the seeds 

The energy momentum tensor is given by the variation of the 
action with respect to the metric. It is 

T,, = d,(3 d,(5 - ^g,, {dxp ■ d'p) . (115) 

The required source functions can be directly derived from this 
expression using (^) to (^). For the scalar sources we use 



$. = ^(/p + 3^/„) ^s = -^s-2U with (116) 



fp = l(P' + {^P?) (117) 



jhj -r — — 

/. = -fe?('3-A.) (118) 

/» = -2^(/9,.-/9,.-3'5..(V/3)^), (119) 

where ^ denotes the Fourier transform (note that (5(5 ^ (5'^ !!). 
Vector sources are determined by 

wf^ =f^,-^^k'M5^i , (120) 

and it is sufficient to calculate the correlation function of one of 
them, e.g. wi, as the transversal character of w imposes 

{wf\k,T)wf\k,T')) = {hkj - k^6,j)V^W{kT,kT'). 
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For the same symmetry reasons the tensor type correlators are 
also determined by one function, F5, alone, see the equation ( |it)i| ) 
in Chapter 0. Therefore we can again pick one special index selec- 
tion to determine F5. Since F5 does not depend on the direction 
of fc, we can choose the special coordinate system ki = k2 = 0, 
leading to F5 = VV? {Ti2{r)Ti2{r')) . 



4.1.3 The large N limit 



Before discussing numerical simulations of global defect, we study 
the limit where the number of components of the scalar field 
becomes very large |55|.. As we shall see, in this limit, the single 
non-linear term in the cr-model can be replaced by its average, 
and the equation of motion becomes linear and can be solved. 
This solution has been found in Ref. [|i4(j|] . Using the equation of 
motion (Ol) in a FLRW metric, = l/a^{(3 - 2a/a(5 - A/3), 
we find 



a 



a ^ 



In a first step we impose scaling on the trace of the energy mo- 
mentum tensor by setting = T{x)/(ar)'^ where T{x) is now 
dimensionless. In the large- approximation, the fluctuations of 
quadratic quantities, like the energy momentum tensor, are of 
order so we neglect them for the field evolution, and we 

replace T{x) with its average T. To solve the resulting equation, 
we change into Fourier space and replace d/ahy a/r, which is ex- 
act in perfectly matter [a = 2) and radiation [a = 1) dominated 
universes, and an acceptable approximation otherwise. This leads 
to 

P + 2^P+ (^e-^^/3 = 0. (122) 

This equation is solved by 

(3{k, t) = r'/'-^iMk) Ukr) + h{k) J-Akr))^ (123) 
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where Jy denotes the Bessel function of order v and v = T -\- 
(112 -af. 



The functions fi{k) are random variables, and we can take them 
to be Gaussian distributed and uncorrelated at all points, 

s:lm 

{fmfri^)) = c\kmk - q)5i^^ . 

We discard the solution with negative z^, since it diverges for 
r — ^ 0. Furthermore, we choose the solution starting as white 
noise as function of fc, leading to n = —2u. Enforcing the a 
model condition that the field cannot leave the vacuum manifold 
((/?(a?,r)^) = 1) fixes C and requires v = a -\- 1. The abso- 
lute normalisation is actually not important for the calculation 
of CMB anisotropics, since the fluctuations will be normalised 
to the COBE data points. In this case C merely determines the 
required energy scale of symmetry breaking. 

Using x(a;) = J„{x)/x'^ as well as (fi{x) = |x(ic) — Jiy+i{x) / x'^~^ 
we can write the solution as 

r) = VAT'/\ikT)Mk) (125) 

and its time derivative as 

r) = VAr'/\{kr)(3Uk) ■ (126) 

It is now easy to derive the energy momentum tensor of the seeds 
using the expressions of the last section. As a worked out example, 
we take a closer look at fp. The equal time correlator (ETC) is 

(/,(fc,r)/;(fc',r)) 

= J d\d^p{(piqT)(p{\k - q\T)(p{pT)(p{\ - k' - p\t) 

-r'^qik - q)p{k' + p)x(gr)x(|fe - q\r)x{pr)x{\k' + p|r) 
+rV(fc' + P)^{qr)^i.\k - q\r)xipr)xi\k' + p\r) 
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-r'^qik - q)(p{pT)(p{\k' + p\T)x{qr)x{\k - g|r)} 
An(Q)An(fc - q)PUp)PU-k' - p)) . (127) 



The expectation value of the initial fields is given by Eq. (|i24| ) and 
the requirement that /3in be a Gaussian random variable implies 



An(g)An(fc - Q)An(p)An(-fc' " P 

C^6{k - k' 



[5{p + q) + 5{p-{q-k))]. (128) 



This allows us to perform the integral over (fip. We introduce the 
dimensionless variables x = qr and y = kr. To simplify the nota- 
tion, we replace all occurrences of an expression like a{x)b{\y—x\) 
by {ab). 



\fp\)iy^r) = ^ J d\ {{^^f+[x{y-x)nxX? 

-2[x{y-x)]{^^){xx)] ^ 

= / dxd^x^ [{^iff + [xy^-x^f{xxf- 

2[xy^x-x\^^){xx)}■ (129) 

In the last equation we performed the integration over one angu- 
lar variable and introduced /i = [x ■ y). 

In this way all required unequal time correlators can be derived. 
They are shown in Figs, g to |7|. A more explicit treatment can be 
found e.g. in [R? . 



4-1-4 Numerical simulation of global texture 

Fields with a finite number of components cannot be treated 
analytically. The unequal time correlators have to be calculated 
numerically on a grid. A useful approach for global field simula- 
tion is to minimize the discretized action |ll8|j . There one does 
not solve the equation of motion directly, but use a discretized 
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version of the action 

S = J (fxa\T) ^d,/3 -9^/3 + ^ - 1) , (130) 

where A is a Lagrange multiplier which fixes the field to the vac- 
uum manifold (this corresponds to an infinite Higgs mass). Tests 
have shown that this formalism agrees well with the complemen- 
tary approach of using the equation of motion of a scalar field 
with Mexican hat potential and setting the inverse mass of the 
particle to the smallest scale that can be resolved in the simu- 
lation (typically of the order of 10"^^ GeV), but tends to give 
better energy momentum conservation. 

As we cannot trace the field evolution from the unbroken phase 
through the phase transition due to the limited dynamical range, 
we choose initially a random field at a comoving time r = 2Ax. 
Different grid points are uncorrelated at all earlier times . 



The use of finite differences in the discretized action as well as 
in the calculation of the energy momentum tensor introduce im- 
mediately strong correlations between neighboring grid points. 
This problem manifests itself in an initial phase of non-scaling 
behaviour, the length of which varies between lOAx and 20Aa;, 
depending on the variable considered. It is very important to use 
results from the scaling regime only (cf. Fig. [I|). 

In order to reduce the time necessary to reach scaling and to 
improve the overall accuracy, one has to choose the finite differ- 
ences in an optimal way. One possibility is to calculate all values 
in the center of each cubic cell defined by the lattice. The addi- 
tional smoothing introduced by this improves energy-momentum 
conservation by several percent. 

To calculate unequal time correlators (UTC), the values of the 
observables under consideration are saved once scaling is reached 
at time Tc and then correlated at all following time steps. While 
there is some danger of contaminating the equal time correlator 
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Fig. 1. The ETCs Cu{z, 1) = (|«'|^)(A;r) (panel a) and C22iz, 1) = (|^'P)(A:r) (panel 
b) are shown for different times. In grid units the times are r = 4 (dashed), r = 8 
(dotted), T = 12 (long dashed), r = 16, 20 (dash dotted, long dash dotted) and 
r = 24 (solid). Clearly C22 scales much sooner than Cu. To safely arrive in the 
scaling regime one has to wait until r ~ 16 and Cij{kT = 0) is best determined at 
r > 20 but kr < 1. 

(ETC), which contributes most strongly to the Q's, with non- 
scahng sources, this method ensures that the constant for kr 
is determined with maximal precision for the ETCs. This is 
very important as the constants Cij(0, 1) fix the relative size of 
scalar, vector and tensor contributions of the Sachs- Wolfe part 
and severely influence the resulting C^'s. In contrast, the CMB 
spectrum seems quite stable under small variations of the shape 
of the UTCs. 

The resulting UTCs are obtained numerically as functions of the 
variables /c, r and Tc with t > Tc and Tc fixed. They are linearly 
interpolated to the required range. One then constructs a her- 
mitian matrix in kr and /cr', with the values of kr chosen on a 
linear scale to maximize the information content, < kr < x^mx- 
The choice of a linear scale ensures good convergence of the sum 
of the eigenvectors after diagonalization (see Fig. g), but still re- 
tains enough data points in the critical region, 0{x) = 1, where 
the correlators start to decay. In practice one chooses as the end- 
point x-aiax of the range sampled by the simulation the value at 
which the correlator decays by about two orders of magnitude, 
typically x^ax ~ 40. The eigenvectors that are fed into the Boltz- 
mann code are then interpolated using cubic splines with the 
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scale [fraction of total size] 



Fig. 2. Energy momentum conservation of numerical simulations is shown. The lines 
represent the sum of the terms which has to vanish if energy (solid) respectively 
momentum (dashed) is conserved, divided by the sum of the absolute values of these 
terms. The abscissa indicates the wavelength of the perturbation as fraction of the 



size of the entire grid, (from [47| 



condition v^^^kr) for kr ^ Xj 



There are several methods to test the accuracy of simulations: 
One of them is energy momentum conservation. In Ref [^7| 



found to be better than 10% on all scales larger than about 4 grid 
units, as is shown in Fig. |2|. Another possibility is a comparison 
with the exact spherically symmetric solution in non-expanding 
space . 



The overall shape and amplitude of the unequal time correlators 
are quite similar to those found in the analytic large- A^^ approxi- 
mation [ |145| , P5| , |45|j (see Figs. || to |7|). The main difference of the 
large- A approximation is that there the field evolution, Eq. ( |i2i| ). 



is approximated by a linear equation. The non-linearities in the 
large-A seeds, which are due solely to the energy momentum 
tensor being quadratic in the fields, are much weaker than in the 
texture model where the field evolution itself is non-linear. There- 
fore, decoherence which is a purely non-linear effect, is much 
weaker in the large-A limit. This is actually the main difference 
between the two models as can be seen in Fig. 0. Otherwise the 
similarity of the results obtained by simulating the full non-linear 
problem and by considering the simplified linear limit is quite re- 
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Fig. 3. The Ci power spectrum is shown for the large-A^ hmit (bold hne) and for 
the texture modeh The main difference is clearly that the large- curve shows some 
acoustic oscillations which are nearly entirely washed out in the texture case. 

markable. The entire class of global 0{N) models behaves in this 
way, and potentially many other global defect models as well. 



4-2 Cosmic strings 



So far we discussed only theories with global defects. Yet in mod- 
ern particle physics, local (gauge) symmetries play a much more 
important role than global symmetries. As we have seen in Chap- 
ter 0, strings are the only local defects which scale and which are 
therefore potential candidates to seed structure formation. Most 
research concentrates on local U{1) theories, the best known cos- 
mic strings. But also other models, even with non-abelian strings, 
have been considered (see e.g. [|iat)| , p^] ). 

The main difference for the numerical treatment of local strings 
as compared to global defects is due to the existance of a gauge 
field which compensates the gradient energy of the field. All the 
energy is therefore concentrated in the tiny region of order the 
symmetry breaking scale Im ~ 1/M ~ l/Tc where the scalar 
field leaves the vacuum manifold. Let us estimate how thin cos- 
mic strings really are: COBE normalisation requires the phase 
transition to take place at the GUT scale, Tc ~ lO^^GeV. This 
corresponds to Im ~ l/^c ~ 10"^^ cm. Clearly, a numerical sim- 
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ulation with a grid of not much more than (512)^ cells which 
should simulate the entire Hubble volume, ~ (lO^^cm)^, cannot 
resolve this scale by more than 55 orders of magnitude. There- 
fore, strings are approximated as infinitely thin, and it can be 
shown that they obey to a very good approximation the Nambu- 
Goto action of fundamental string theory [|152|] . Corrections are of 
the order of the string thickness devided by the string curvature 
scale, and therefore irrelevant for cosmology. 

Furthermore, the string network needs to loose energy by grav- 
itational radiation in order to scale. Hence, energy momentum 
conservation cannot be enforced for the defect field alone, lead- 
ing to 14 UTCs insted of only five. As the interactions of cosmic 
strings are not known, one must also make ad hoc assumptions 
concerning the decay products. This choice has considerable im- 
pact on the results [[25| , |i24| . 



We do not describe the numerical simulations to evolve cosmic 
strings. Detailed accounts of this problem can be found in the lit- 
erature [[r52| , |7T| , pXf6| , p^j . Let us, nevertheless, point out the main 



problem. It is very difficult to simulate cosmic strings in expand- 
ing space, due to the large difference between the Hubble scale 
and the scale of small scale structure. Hence, it remains unclear 
up to date, whether string simulation in an expanding universe 
can capture enough of the small scale structure, the tiny wiggles 
and loops which develop due to the string self-interaction, to pro- 
duce meaningful results. On the other hand, string simulations in 
flat space, do not satisfy energy momentum conservation of ex- 
panding space. To adress the small scale structure problem, most 
of the recent results in the literature, actually all except [^,0,|7j], 
use flat space simulations or semi-analytical methods to calculate 
the string UTCs. It is not clear to us which procedure gives the 
best results, but since all the obtained CMB spectra disagree sig- 
niflcantly with observations, this question has somehow lost its 
urgency. 

As initial conflguration of a string simulation, one usually lays 
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down string segments according to the so called Vachaspati-Vilenkin 
algorithm. These are then evolved with the Nambu-Goto equa- 
tion of motion. The physical problem of the 'decay product' of 
cosmic strings is related to the numerical problem of small scale 
structure: Most string codes find that the network develops struc- 
ture (wiggles, tiny loops) on the smallest scales which the simu- 
lation can resolve. The phyical scale of these small wiggles and 
loops is still unknown. It may even be, that the loops become 
smaller and smaller due to self-intersection, until their size is of 
the order of their thickness and they decay into elementary par- 
ticles. This picture, which is in contrast to the decay into gravity 
waves, is avocated in Ref. f|i54J]. There have also been several at- 



tempts to take into account these wiggles in semi-analytic mod- 
els im . 
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5 Result 



5.1 The unequal time correlators 



As explained in Chapter g to compute the observable CMB and 
matter, power spectra we need the unequal time correlators of 
the seed energy momentum tensor. 

More precisely, for the scalar part we need the correlators 

r)$:(fc, t')) = ^^=Cn{z. r) , (131) 

($,(fc, r)^:(fc, t')) = ^^Cu{z. r) , (132) 

(*,(fc, r)*:(fc, t')) = ^^C22{z. r) , (133) 

as well as C2i{z,r) = Cl2{z, 1/^)- The functions Cij are analytic 
in z'^. The pre-factor 1/ {k'^^/rr') comes from the fact that the 
correlation functions (fpfp), k^{fT^fl) and k'^{fvfl) have to be 
analytic and from dimensional considerations (see Ref. [|^]). 



The functions Cij are shown in Fig. Panels (a) are obtained 
from numerical simulations. Panels (b) represent the same corre- 
lators for the large- limit of global 0(A)-models (see [pr^ , |95|j ) . 



In Fig. 1^ we show Cij{z, r = 1), and the 'constant' of the Taylor 
expansion for Cij is given as a function of r, i.e., Cij(0,r). 

Vector perturbations are induced by E*-*^ which is seeded by w^^\ 
Transversality and dimensional arguments require the correlation 
function to be of the form 

{wt'\k,t)wf^\k,T')) = v^(A;% - k,k^)W{z,r) . (134) 
Again, as a consequence of causality, the function W is analytic 
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0.01 



Fig. 4. The two point correlation functions Cii{z,r) = rT'(<I>s(fc, t)^l(k, r')) 
(left), C22{z,r) = k'^VT?{^s{k,T)^*s{k,T')) (center) and 

\Ci2{z,r)\ = k^VTT'\{^s{k,T)'i'l{k,T'))\ (right). Panels (a) represent the re- 
sult from numerical simulations of the texture model; panels (b) show the large- A'^ 
limit. For fixed r the correlator is constant for z < 1 and then decays. Note also 
the symmetry under r ^ 1/r for Cu and C22 which is lost for C12 (from j47||). 



in (see The functions W{z,r), W{z,l) and W{0,r) are 

plotted in Fig. ||. 



Symmetry, transversality and tracelessness, together with sta- 
tistical isotropy require the tensor correlator to be of the form 
(see [|5|]) 



'^ij^lm H" k {Sijkikm + Sifjikikj dukjkm ^imkik 
^jlkikjji Sjm 



5jikikra ^jmkiki) + k ^kikjkikm] ■ (135) 



The functions T{z, r) as well as T{z^ 1) and r(0, r) are shown in 
Fig. 0. 
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Fig. 5. On the left the correlators Cij{z, 1) are shown, while the right figure depicts 
Cij{0,r) with r = t'/t. The solid, dashed and dotted lines represent C22 , Cn 
and IC12I respectively. Panels (a) are obtained from numerical simulations of the 
texture model and panels (b) show the lavge-N limit. A striking difference is that 
the large- value for IC12I is relatively well approximated by the perfectly coherent 
result y^|CiiC22| while the texture curve for IC12I lies nearly a factor 10 lower 
(from 



Expanding the unequal time correlators in terms of eigenfunc- 
tions and eigenvectors as explained in Chapter 0, one finds that 
ae expansion with a linear weight actually converges faster than 
one with a logarithmic weight. This is illustrated in Fig. g below. 

In Fig. we also show the equal time correlators of the energy 
density and pressure, velocity and anisotropic stresses of local 
cosmic strings. It is interesting to note to which extent the spatial 
components of the energy momentum tensor are smaller than the 
energy density. In contrary to global defects, local strings provide 
a nearly non-relativistic source. 



A source is called totally coherent [|iU5| , |^] if the unequal time 



correlation functions can be factorized. This means that only one 
eigenvector is relevant. A simple totally coherent approximation, 
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(a) 




Fig. 6. The vector correlator W{z, r) is shown on the left. The texture simulations, 
panel (a), and the \aige-N limit, panel (b), give very similar results, up to a slight 
difference in amplitude. On the right we see W{z,l) (top) and VF(0, r) (bottom). 
The solid line represents the texture simulations and the dashed line is the large- 
result. The 'wings' visible in the lower texture curve are probably not due to a 
resolution problem but the beginning of oscillations (from p7|). 

which however misses some important characteristics of defect 
models, can be obtained by replacing the correlation matrix by 
the square root of the product of equal time correlators, 

(S,(r)S-(r')> - ±v/(|5,(r)P)(|5,(r')P> . (136) 

This approximation is exact if the source evolution is linear. Then 
the different k modes do not mix and the value of the source term 
at fixed k at time r is given by its value at initial time multiplied 
by some transfer function, <S(fc,r) = T{k^r,Tin)S{k,rin). In this 
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(a) 






Fig. 7. The tensor correlator T{z, r) is shown on the left, the texture simulations 
in panel (a), and the large- limit in panel (b). On the right we see T(z, 1) (top) 
and T(0, r) (bottom). The solid lines represent texture simulations and the dashed 
lines are the large- result, (from 



situation, ( |i5t3|) becomes an equality and the model is perfectly 
coherent. Decoherence is due to the non-linearity of the source 
evolution which induces a 'sweeping' of power from one scale into 
another. Different wave numbers k do not evolve independently. 

It is interesting to note that the perfectly coherent approximation 
(|i5b|) leaves open a choice of sign which has to be positive if i = 
J, but which is undetermined otherwise. According to Schwarz 
inequality the correlator (<Si(r)<S*(r')) is bounded by 



"WmiiWW) < {Si(r)SUT')) < ^{\S,(rW){\Sj(T')P] 



Hence, for scales/ variables for which the Greens function is not 
oscillating (e.g. Sachs Wolfe scales) the full result always lies be- 
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Fig. 8. The sum of the first few eigenfunctions of T{x, x) is shown for two differ- 
ent weight functions, (a) logarithmic, w = 1/x and (b) linear, w = 1. The first 
(long dashed), first and second (short dashed), first ten (dotted) and first thirty 
(sofid) eigenfunctions are summed up. The open circles represent the full corre- 
lation function. Clearly, the eigenfunctions obtained by linear weighting converge 
much faster. Here we only show the equal time diagonal of the correlation matrix, 
but the same behavior is also found in the Ci power spectrum which is sensitive to 
the full correlation matrix (from Re. [p7|] ). 
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Fig. 9. The size of the equal time contribution from the energy density (solid line) 
and the anisotropic stresses (left). Pressure, velocity and vorticity are of about the 
same size as the anisotropic stresses (right) (from Magueijo & Brandenberger [[106|). 



tween the 'anti-coherent' (minus sign) and the coherent result. 
This behavior has been verified numericaUy [j57 . 



The first evidence that acoustic peaks are suppressed in defect 
models has been obtained in the perfectly coherent approxima- 
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tion in Ref. [j?^. In Fig. [TD|we show the contributions to the C^'s 
from more and more eigenvectors. A perfectly coherent model has 
only one non-zero eigenvalue. 




Fig. 10. The scalar, vector and tensor contributions for the texture model of struc- 
ture formation are shown. The dashed lines show the contributions from single 
eigenfunctions while the solid line represents the sum. Note that the single contri- 
butions to the scalar and tensor spectrum do show oscillations which are however 
washed out in the sum. (Vector perturbations do not obey a wave equation and thus 
do not show oscillations.) 

A comparison of the full result with the totally coherent approx- 
imation is presented in Fig. [TT]. There one sees that decoherence 
does smear out the oscillations present in the fully coherent ap- 
proximation, and does somewhat damp the amplitude. Decoher- 
ence thus prevents the appearance of a series of acoustic peaks. 
The absence of power on the angular scale of the peaks, how- 
ever, is not a consequence of decoherence. Is mainly due to the 
anisotropic stresses of the source, which lead to perturbations in 
the geometry inducing large scale Q's (Sachs Wolfe), but not to 
density fluctuations, which are responsible for the acoustic peaks. 
Large anisotropic stresses are also at the origin of vector and ten- 
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Fig. 11. The power spectrum for the texture scenario is shown in the perfectly 
coherent approximation (top panel) and in the full eigenfunction expansion. Even 
in the coherent approximation, the acoustic peaks are not higher than the Sachs 
Wolfe plateau. Decoherence just washes out the structure but does not significantly 
damp the peaks. 

sor fluctuations which contribute more than 50% on large scales. 
These results were first computed fully in Ref. [|iiY|] and followed 
by a parameter study in Ref Thy are in agreement with 
Ref. but disagree strongly with Ref. []2§ . 



In the real universe, perfect scaling of the seed correlation func- 
tions is broken by the radiation-matter transition, which takes 
place at the time of matter and radiation equality. Teg ~ 20/i^^Q^^/^Mpc. 
The time Tgg is an additional scale which enters the problem and 
influences the seed correlators. Only in a purely radiation or mat- 
ter dominated universe are the correlators strictly scale invariant. 
This means actually that the k dependence of the correlators C, 
W and T cannot really be cast into a dependence on x and x', but 
that these functions depend on r, r' and /c in a more complicated 
way. In principle, one thus has to calculate and diagonalize the 
seed correlators for each wave number k separately and the huge 
gain of dynamical range is lost as soon as scaling is lost. 

In the actual case at hand, however, the deviation from scaling 
is weak, and most of the scales of interest for structure forma- 
tion (especially those which are still in the linear regime) enter 
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the horizon only in the matter dominated era. The behavior of 
the correlators in the radiation dominated era is of minor impor- 
tance. To solve the problem, one can calculate the eigenvalues 
and eigenfunctions twice, in a pure radiation and in a pure mat- 
ter universe and interpolate the source term from the radiation 
to the matter epoch. Denoting by Xm^Vm and Xr^Vr a given pair 
of eigenvalue and eigenvector in a matter and radiation universe 
respectively, we choose as our deterministic source function 



v{r)=y{r)\/XrVr{kr) + (1 - y{r)) \/ X„,Vm{kr) (137) 
with, e.g., 

■J- 

y{r)= ^'^ or y{r) = exp{-r/req) , (138) 

r + Teq 

or some other suitable interpolation function. The effect of the 
radiation dominated early state of the universe is relatively unim- 
portant for the scales of interest for CMB anisotropics and linear 
gravitational clustering The difference between the pure matter 
era result and the interpolation is very small This seems 



to be quite different for cosmic strings where the fluctuations in 
the radiation era are about twice as large as those in the matter 
era [|132|j . The radiation dominated era has, however, little effect 
on the key feature of CMB anisotropics from topological defects; 
namely the absence of acoustic peaks. 

In models with cosmological constant, there is actually a second 
break of scale invariance at the matter-A transition. There one 
can proceed in the same way as outlined above. Since defects 
cease to scale and disappear rapidly in an exponentially expand- 
ing universe, the eigenvalues for the A dominated universe all 
vanish. 



In Fig. |r2| we show the CMB power spectrum from cosmic strings. 



Predictions from different research groups working on the sub- 
ject [|5|,[7|, [251 , pi]51 , |11^4| , [r7Tl] agree only partially. Contrary to global 



defects models, there seems to be a broad peak at rather high 
£ ~ 400 to 500. The height of the peak is very much a matter of 
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debate. In some work it is completely absent [j^, while in other it 
is quite high [|2^ , |il^ij| . We can maybe understand this difference 
between global and local defects as being due to the larger dif- 
ference between the source functions in the radiation and matter 
dominated era: For the global 0{N) defects, the ratio between 
the sources in the two epochs is about 1.2, while for strings it is 
rather in the vicinity of 4 [ |lb(J|j . This difference can explain the 
more prominent peak at high i for cosmic strings. An additional 
intriguing difference to global defect is the strong domination of 
the energy density over the rest of the energy momentum tensor 
(see Fig. ^. The dependence of the spectrum on the equation of 
state of the decay product, px = wxPx is remarkable. 
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Fig. 12. The angular CMB power spectrum of cosmic strings, from a simulation 



by Contaldi, Hindmarsh and Magueijo |26|. The figure shows the result for strings 
decaying into a fluid obeying the equation of state p = wp. with w = 1/3, 0.1, 0.01 
respectively. The old data overlaid to the graph is to be ignored. 



A feature which the two classes of defects have in common is 
decoherence which smears out all secondary peaks. Furthermore, 
cosmic strings seem to suffer as well from insufficient power on 
large scales in the dark matter. The clear sequence of acous- 
tic peaks recently observed [|ilU| . |BB| , P7|] rules out any seeds with 
strongly non-linear evolution, such as 0{N) (with low A^^), both 
for global and local gauge theories, as well as any other non-linear 
mechanism for seeding cosmic perturbations! 
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5.2 non-Gaussianity 



An interesting difference between structure formation with topo- 
logical defects and inflationary models is also that the latter have 
generically Gaussian perturbations while the former don't. Even 
if the defect energy momentum tensor would be Gaussian ini- 
tially, non-linear evolution induces non-Gaussianities. (Recently, 
however, a model leading to Gaussian fluctuations from cosmic 
defects has been investigated 1^.) 

Even if the fields themselves would be Gaussian (which they are 
not), the energy momentum tensor, which is quadratic in the 
fields, would obey a x^-distribution. But the cr-model condition 
Eili(A)^ = 1 cannot be satisfied if the fields /3i are Gaussian 
random variables. In the large model, the 0(A)-model in the 
limit A ^ oo, the energy momentum tensor is a infinite sum 
of variables d^Pid,y(3i which all obey the same distribution. This 
large A model is nearly Gaussian as a consequence of the central 
limit theorem. The fact that the variables Pi have to satisfay 
the normalization condition however implies that the variables 
d^PidjyPi are not statistically independent. 

Clearly, 0(A)-models with A > 4 components do not lead to 
topological defects in 3-dimensional space, but in the large A- 
limit the equations of motion can be solved analytically and, as 
we have already seen, the model is very useful to study certain 
features of 0(A) defects. The fact that non-Gaussianity becomes 
weaker as A becomes larger, indicates that cosmic strings are 
probably the most non-Gaussian defects and textures are the 
least non- Gaussian. 

For defects to be non-Gaussian, means that their energy momen- 
tum tensor does not obey Gaussian statistics. Therefore certain 
reduced higher oder moments do not vanish. It is not evident 
how to find the best observable to locate the non-Gaussianity. 
Some suggestions for variables which might be useful in the case 
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of defects have been proposed in Refs. [j54|,B9|. Furthermore, a 



given observational variable like, e.g. , the integrated Sachs- Wolfe 
effect, may be the sum of many non-Gaussian but equally dis- 
tributed contributions and hence be very closely Gaussian due to 
the central limit theorem. 

Unfortunately, it is not known how strong non Gaussian features 
are in the CMB or in the dark matter distribution of topolog- 
ical defect models. It is also not clear on which scales they are 
strongest. Due to the arguments indicated above, we expect cos- 
mic strings to be most non-Gaussian. A well distinguished non- 
Gaussian feature in CMB anisotropics from cosmic strings is the 
Kaiser-Stebbins effect [^. This is a disconinuity in the CMB 
temperature due to a moving string between the observer and 
the CMB. Even though this effect is easily obtained analytically 
for a straight cosmic string, it is difficult to get a handle on it in a 
string network of many bent and twisted strings which also con- 
tain small scale structure. One expects the Kaiser-Stebbins effect 
to be relevant for the anisotropics on about arc-minute scale. 

Old partial results on non-Gaussianity from numerical simula- 
tions (see e.g. . [|ii8| , |5T|j ) are probably strongly affected by finite 



size effects and not very reliable. On the other hand, with the 
more successful method to compute power spectra just by de- 
termining the unequal-time two point distribution of the defects, 
one loses all information about higher order correlations and thus 
about non-Gaussianity. Since the known defects are such a bad fit 
to the present CMB data, nobody has been sufficiently motivated 
to study and solve the difficult problem of the non-Gaussianity 
which defects may induce in the CMB anisotropics or in the 
matter distribution. A semi-analytic study of the bi-spectrum, 
the third moment of CMB anisotropics for cosmic strings can be 
found in Ref. [j5^ . 



We know that also non-linear Newtonian clustering induces non- 
Gaussianities. The reduced n-point functions in the matter dis- 
tribution due to non- linear clustering scale like D^^^^, where D 
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Fig. 13. The power spectrum is shown for different values of cosmological pa- 
rameters. In the right panel we choose f^A = 0, ^cdm = 0.95, Q.^ = 0.05 and vary h. 
In the left panel we fix /i = 0.5, Vlb = 0.05 and vary Qa- We only consider spatially 
flat universes, f^o = 1- 



is the matter density perturbation amplitude (see e.g. [)5^). In 



contrary, a non-Gaussianity already present at the linear level 
simply scales like in the n-point function Therefore, ini- 



tial non-Gaussianities, hke those of defects are best detected on 
large scales, where non-linearities are negligible since D <C 1. 
On small scales, where non-linearities are important, the non- 
Gaussianity from Newtonian clustering is stronger than the one 
form the initial conditions. Therefore, the fact that non-linearities 
in the observed galaxy distribution closely follow the behavior ex- 
pected from Newtonian clustering, does not constrain topological 
defects very strongly. 



5. 3 CMB anisotropies and polarization 



The C/s for the global texture model are shown in Fig. [TT], bot- 
tom panel. 

Vector and tensor modes are of the same order as the scalar com- 
ponent at COBE-scales. For the 'standard' texture model one 
has d? : C\V : c\V - 0.9 : 1.0 : 0.3. These results were found 
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numerically and by analytical arguments in Refs. |T17|,5],E],ff7| and 



^j . Due to tensor and vector contributions, even assuming per- 



fect coherence (see Fig. [TT], top panel), the total power spectrum 
does not increase from large to small scales. Decoherence leads to 
smoothing of oscillations in the power spectrum at small scales 
and the final spectrum has a smooth shape with a broad, low 
isocurvature 'hump' at ^ ~ 100 and a small residual of the first 
acoustic peak at £ ~ 350. There is no structure of peaks at small 
scales. The power spectrum is well fitted by the following fourth- 
order polynomial in x = logi: 

llOCio 



The effect of decoherence is less important for the large- model, 
where oscillations and peaks are still visible (see Fig 0, bottom 
panel). As argued before, this is due to the fact that the non- 
linearity of the Idiige-N limit is only in the quadratic energy 
momentum tensor. Since decoherence is inherently due to non- 
linearities, we expect it to be stronger for lower values of A^. 

In Fig. [T3| we plot the global texture Cg power spectrum for dif- 
ferent choices of cosmological parameters. The variation of pa- 
rameters leads to similar effects like in the inflationary case, but 
with smaller amplitude. At small scales {£ > 200), the Cgs tend 
to decrease with increasing Hq and they increase when a cos- 
mological constant Qj^ = 1 — is introduced. Nonetheless, the 
amplitude of the anisotropy power spectrum at high is remains 
in all cases on the same level like the one at low is, without 
showing the substantial peak found in inflationary models and in 
the data. The absence of acoustic peaks is a stable prediction of 
global 0{N) models. The models are normalized to the full CMB 
data set, which leads to slightly larger values of the normahzation 
parameter than pure COBE normalization. An avarage value for 
different choices of cosmological parameters is e = AnGrf ~ 1.6 

We compare the texture results with the three best current ex- 
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perimental data sets from Boomerang-98 [|iiU|] , Maxima and 



DASI We take in to account the calibration uncertainty of 



10% for Boomerang and of 5% for Maxima and DASI. We also 
include the COBE dataset using Lloyd Knox's RADPack pack- 
ages [|i(3(3 . 



No matter how precisely an experiment will measure the CMB 
sky, the result will be always affected by the intrinsic statistics 
of the perturbations, i.e. the cosmic variance. In the recovery of 
the power spectra of the above experiments, the temperature dis- 
tribution is often assumed to be Gaussian. Deviation from Gaus- 
sianity could lead to an enhancement of the cosmic variance, 
which can be as large as a factor of 7 (see [|iu;j|] ). Inevitably 



therefore, the quoted error-bars in the CMB power spectrum 
could be in principle underestimated by the assumption of the 
Gaussian statistic. This effect is, however only relevant for rela- 
tively low is. Furthermore, preliminarly statistical analysis of the 
Maxima map [|i61|j have found no-evidence for non-Gaussianities. 
Keeping this caveat in mind, but missing a more precise alterna- 
tive, we indicate the minimal, Gaussian error calculated accord- 
ing to the published data. The numerical seeds taken from 



are assumed to be about 10% accurate. In Fig. [T^ we plot the 
Boomerang and DASI data together with the theoretical predic- 
tions for a texture model with h = 0.65, Ha = 0.625, Qcdm = 0.315 
and fib = 0.06, and the corresponding inflationary model with 
Us = 0.98 which gives a best fit to the data. All of the exper- 
iments detect a very clear first peak at £ « 200 which is in- 
compatible with all 0{N) models [p|, pT7| , rr2g| , |47|] as well as with 
the cosmic string results [^,|T2T],|T2^ In Table we report the 
values from a comparison of the CDM-texture model with 
each experiment, separately. As we can see, while in reasonable 
agreement with the large angular scale data from the COBE ex- 
periment, the new data at intermediate and small angular scales 
from BOOMERanG, DASI and MAXIMA rule out the model 
with high significance. 
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Table 1 

Topological defects versus CMB data. While in reasonable agreement with the 
COBE data, the power spectrum inferred from textures is in strong disagree- 
ment with the present intermediate and small angular scale observations. The 
BOOMEranG data alone, in particular, rule out the model at extremely high sig- 
nificance. 

In addition, the data shows indications of a series of peaks point- 
ing to coherent acoustic oscillations. If the evidence for this peak 
sequence strenghtens once new data becomes available, it will 
provide strong evidence against non-linear mechanisms for the 
generation of perturbations. Furthermore, it will enable us to 
place stringent limits on the relative amounts of isocurvature and 
adiabatic contributions to the perturbations. We present a more 
general discussion in the next chapter. 

Finally, we compare E- and 5-type polarisation induced by de- 
fects with the result of an inflationary model with equal ampli- 
tude of scalar and tensor temperature fluctuations in Fig. [T5] . Due 
to the amount of vector perturbations present, it is not surpris- 
ing that 5-type polarization has significantly higher amplitude 
for defects. Furthermore, the isocurvature shift of the 'acoustic 
peaks' and decoherence are also visible in the polarisation signal. 



5.4 Matter power spectra 



The first serious attempts to compute the matter distribution in 
scenarios with topological defects were done for cosmic strings [pXfS 
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Fig. 14. The Ci spectrum obtained in the standard texture model is compared with 
the recent BOOMERANG and DASI data. A standard inflationary model, which 
gives a much better fit to the data, is also plotted. 
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Fig. 15. Power spectra of temperature (T), electric type polarization (E) and mag- 
netic type polarization (B) for global strings, monopoles, textures and an 0{N) 
scalar field with > 4 are shown. For comparison, also the corresponding spectra 
in a standard CDM model with T/S = 1 (which maximizes the B component) is 
plotted. All defect models predict a much larger component of B polarization on 
small angular scales (from Seljak, Pen and Turok |128| ). 
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But due to the problems inherent in cosmic string simulations, 
the accuracy of these results may well be less that a factor of two 
or three. 



The situation is different for global defects. There, the first ac- 
curate matter spectra have been computed in Ref. [ |117|j . This 
work, which underestimated the vector and tensor contribution 
to CMB anisotropics, still had the wrong normalization. This was 
corrected later in [[5T| , |iiY|p7| . For a COBE mormalized global 
texture model, the total mass fluctuation aji in a ball of ra- 
dius R = 8/i"^Mpc, is about ag = (0.44 ± 0.07)/i (the error 
coming from the CMB normalization) in a flat universe with- 
out cosmological constant (see Refs. [|iiY| , |47| . From the observed 



cluster abundance, one infers dg = (0.50 ± 0.04)^ and 



C^8 



o.59i:S 



P^] . These results, which are obtained with the 
Press- Schechter formula, assume Gaussian statistics. We thus 
have to take them with a grain of salt, since we do not know 
how non-Gaussian fluctuations on cluster scales are in the tex- 
ture model. According to Ref. [j57j, the Hubble constant hes in 
the interval h ~ 0.73 ±0.06 ±0.08. Hence, in a flat CDM cosmol- 
ogy, taking into account the uncertainty of the Hubble constant, 
the texture scenario predicts a reasonably consistent value of ag. 



However, as discussed in Refs. [j2[, [|iiY|] and unbiased global 
texture models are unable to reproduce the power of galaxy clus- 
tering at very large scales, ~ 20h~^ Mpc. In order to quantify 
this discrepancy we compare the prediction of the linear matter 
power spectrum with the decorrelated linear power spectrum of 
the PSCz survey from [j57j in Fig. \n\. Clearly, the textures model 
without a cosmological constant provides a bad fit to the large 
scale data. Including a cosmological constant improves the agree- 
ment between the shape of the theoretical spectrum and the data. 
However, the required bias is too high, 6 ~ 6 and the model fails 
to match the observed value of erg. We further test this discrep- 
ancy by comparing the theoretical predictions with the results 
from a wide number of infrared (Refs. [[56|j , [|i42|j ) and optically- 
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selected (Refs. [|r02|) galaxy redshift surveys, and with the 



real-space power spectrum inferred from the APM photometric 
sample (Ref. [|TD|)- AH the results are in rough agreement with 
the PSCz comparison: For h = 0.5 and Qa = 0.8 the shape of the 
texture power spectrum fits the data well, but the bias is very 
high. The only exception is the APM data, the shape of which 
cannot be fitted with the texture spectrum. In Table |2| we report 
for each survey the best cosmological parameter (/i, Qa) within 
the limited ones which have been analysed, and the best value of 
the bias parameter obtained by x^-minimization. We also indi- 
cate the value of (not divided by the number of data points) 
and the number of data points. The bias parameter strongly de- 
pends on the data considered. This is not surprising, since also 
the catalogs are biased relative to each other. A more detailed 
discussion of this analysis is found in Ref. . 

The power spectra for the large-A^^ limit and for the coherent 
approximation are typically a factor 2 to 3 higher (see Fig. [T^), 
and the biasing problem is alleviated for these cases. For = 
we find ag = 0.57/i for the large- N limit and ag = 0.94/i for 
the coherent approximation. This is no surprise since only one 
source function, the analog of the Newtonian potential, seeds 
dark matter fluctuations and thus coherence always enhances the 
unequal time correlator. The dark matter Greens function is not 
oscillating, so this enhancement translates directly into the power 
spectrum. 

Models which are anti-coherent in the sense defined in Section IID 
reduce power on Sachs- Wolfe scales and enhance the power in 
the dark matter. Anti-coherent scaling seeds are thus the most 
promising candidates which may cure some of the problems of 
global 0{N) models. 

The analysis described here does not take into account the effects 
of non-linearities and redshift distortions. Redshift distortions in 
the texture case should be less important than in the inflationary 
case since the peculiar velocities are rather low (see next para- 
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Fig. 16. The dark matter power spectrum for the texture model (soUd hne) is 
compared with the coherent approximation (short dashed) and the large- hmit 
(long dashed). The spectra are COBE normahzed and the cosmological parameters 
are Q\ = , h = 0.5. 

graph). Non-linearities typically set in at /c > 0.5/iMpc^^ and 
should not have a big effect on our main conclusions which come 
from much larger scales. Inclusion of these corrections will result 
in more small-scale power and in a broadening of the spectra, 
which even enhances the conflict between models and data. Fur- 
thermore, variations of other cosmological parameters, like the 
addition of massive neutrinos, hot dark matter, which are not 
considered here, will result in a change of the spectrum on small 
scales but will not resolve the discrepancy at large scales. 



5.5 Bulk velocities 



To confirm the missing power on 20 to 100/i^^Mpc, we also con- 
sider the velocity power spectrum which is not plagued by the bi- 
asing problems of the galaxy distribution. The assumption that 
galaxies are fair tracers of the velocity field seems to us much 
better justified, than to assume that they are fair tracers of the 
mass density. We therefore test the global texture model against 
bulk velocity data (see also [j?7]). We use the data by Ref. 
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Fig. 17. The matter Power spectrum: comparison between the PSCz data and 
theory. 
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Table 2 

Analysis of the matter power spectrum. In the first column the catalog is indicated. 
Cols. 2 and 3 specify the model parameters. In cols. 4 and 5 we give the bias 
parameter inferred by minimization as well as the value of x^. Col. 6 shows the 
number of 'independent' data points assumed in the analysis. 

which gives the bulk flow 

^1{R) = / P{k)W{kR)dk , (140) 



in spheres of radii = 10 to 60/i^^Mpc. These data are derived 
after reconstructing the 3— dimensional velocity field with the 
POTENT method (see and references therein). 
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Table 3 

Bulk velocities: Observational data from |36| and theoretical predictions. A^, es- 
timates the observational uncertainty. The uncertainties on the theoretical predic- 
tions are around ~ 30%. The models JIa = with h = 0.5 and h = 1 as well as 
JIa = 0.8, h = 0.5 are given. 

As shown in Table 0, the COBE normahzed texture model pre- 
dicts too low velocities on large scales when compared with PO- 
TENT results. More recent measurements of the bulk flow lead 
to somewhat lower estimates like (7y{R) ~ (230 ± 90) oX R = 
60/i~^Mpc (see Ref. [^), but a discrepancy of about a factor 
of 2 in the best case remains. Including a cosmological constant 
helps at large scales, but decreases the velocities on small scales. 
Similar results are obtained in models with cosmic strings. There, 
however, details depend on assumption about the fluid into which 
string loops decay. For loops decaying into radiation, models with 
cosmic strings have even less power on 50/i^^Mpc. 



This comparison of models with data shows that U{1) cosmic 
strings and global 0{N) models are clearly ruled out by present 
large scale structure data. In the next chapter we shall study 
whether generalizations or mixed models may survive. 
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6 Generalizations 



So far we have studied some specific models of topological defects 
and we have found that they cannot reproduce the observed large 
scale structure of the universe. Even if they basically lead to a 
Harrison-Zel'dovich spectrum on large scales, the 'details' just 
do not fit. The isocurvature nature of the induced perturbations 
together with the high amplitude tensor and especially vector 
contributions, lead to relatively high amplitude of fluctuations 
on very large scales. Therefore, the COBE-normalized acoustic 
'peaks' and the dark matter power spectrum typically are too 
low. In addition, decoherence smears out the anyway too low 
secondary peaks in the CMB anisotropy spectrum. 

In this chapter we first study the question whether this feature is 
generic or just happens to occur in the specific models which we 
have analyzed in the preceding chapters. Then we want to investi- 
gate whether mixed models with inflationary perturbations and a 
defect component might fit the data better than pure inflationary 
perturbations. 

6. 1 Generic properties of the unequal time correlation functions 

We start by deriving some generic properties of power spectra 
induced by topological defects. As we have seen in Chapter 2, 
the resulting perturbation power spectra are fully determined by 
the un-equal time correlators of the defect energy momentum 
tensor. 

The main properties which we shall use to describe the defect cor- 
relators are causality, statistical homogeneity and isotropy and 
scaling. Since the nature of the defects resulting as a 'topologi- 
cal relicts' from a phase transition does not enter directly in this 
analysis, we speak of 'scaling causal seeds', seeds denoting any 
non-uniformly distributed form of energy which provide a pertur- 
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bation to the homogeneous background fluid. In first order per- 
turbation theory they evolve according to the unperturbed (in 
general non-linear) equations of motion. For simplicity, we as- 
sume the seeds to be coupled to the cosmic fluid only via gravity. 
A counter example to this are the U{1) cosmic strings discussed 
before. As we have seen, the string CMB anisotropy power spec- 
trum, especially the height of the acoustic peak, depends very 
sensitively on the details of the coupling of string seeds to mat- 
ter [|26| , |i:^4|] . For uncoupled seeds the energy momentum tensor is 
covariantly conserved. To determine power spectra or other ex- 
pectation values which are quadratic in the cosmic perturbations, 
we just need to know the unequal time correlation functions of 
the seed energy momentum tensor [|451 , ^] , 



(e^,(fc, r)e:^(fc', /)) = M%,,,{k, r, r')6{k - k') , (141) 

where M is a typical energy scale of the seeds {e.g. the sym- 
metry breaking scale for topological defects), which determines 
the overall perturbation amplitude. Here we define Fourier trans- 
forms with the normalization 

/(k) = ^ / d'xf{^) exp(zkx) ; (142) 



/(x) = ^ / Skf{\^) exp(-^kx) . (143) 



(27r) 

Seeds are causal., if 



C^^^p(x, r, t') = -^{e^^iy, r)eip{y + x, t')) 



vanishes for |x| > r + r'; and they are scaling., if C depends on no 
other dimensional parameter than fc, r and t' . From scaling we 
conclude that for purely dimensional reasons, we can write the 
correlations functions in the form 

C^,Ap(k, r, t') = ^F^,Ap(v^ • k, r'/r) , (144) 
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where F^^xp is a dimensionless function of the four variables Zi = 
V T'rki and r = r'/r, which is analytic in zi. 



We also require the seed to decay inside the horizon, which implies 
lim C^.Ap(k, r, r') = lim C^.xpiK r, r') = . (145) 

fcr— >oo fcr'^oo 

Furthermore, since the seeds interact with the cosmic fluid only 
gravitationally, O satisfies covariant energy momentum conserva- 
tion, 

ef," = . (146) 

With the help of these four equations, we can, for example, ex- 
press the temporal components, Gq^ in terms of the spatial ones, 
Qij. The seed correlations are thus fully determined by the spa- 
tial correlation functions Cijim. Statistical isotropy, scaling and 
symmetry in i,j and /,m as well as under the transformation 
k;T ^ I, rri] —k] t' require the following form for the spatial 
correlation functions: 

CijlmiK T, t') = 

—=[z^ZjZiZmFi{z^, r) + 

+ ZjZidim + ZjZm5il)F2{z^ , v) + 

ZiZj5imF^{z'^ , r)/r + ziZmSijFz{z^ , l/r)r + 

+Sij6imF^{z^ r) + {6uSjm + SimSji)F5iz^, r)] , (147) 

where the functions Fa are analytic in = rr'k'^, and for a 7^ 3 
they are invariant under the transformation r — > 1/r, Fa{z'^, r) = 
Fa{z'^ ,1/r). The positivity of the power spectra Cijij(k,r,T) = 
> leads to a series of positivity conditions for the func- 
tions Fa'. 

0<F,{z^l), (148) 
0<F4(z^l) + 2^5(^2^ 1) , (149) 
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Q<z'F2{z\l) + Fr,{z\l) , (150) 

0<^^Fi(22,l) + 4^2^2(22,1) + 3^5(^2,1) , (151) 
< z^Fi{z\ 1) + 2z\Fi{z\ 1) + 2^2(^2, 1)) 

+^4(^2,1) + 2^5(^2,1) . (152) 

Since Cjj/^ is the Fourier transform of a real function, 

Cij7m(k, r, t'Y = C^j7m(-k, r, t') . (153) 

Hence the ansatz ( |14V|) imphes that the functions Fa{z'^,r) are 
real. Furthermore, decay inside the horizon (condition ( |145|) ) 
yields 

lim Fa{z\r) = lim Fa{z\r) = . (154) 

z'^r—^oo z^ jr—tco 

In addition, analyticity implies that the functions F^ do not di- 
verge in the limit z ^ 0, thus 

X\mFa{z^,r) = A^ir) 

with 

Aair) = Aa{\lr) for aU a 7^ 3 . 

As an example, we show some of these functions in the large 
limit of global scalar field seeds (see Chapter g). In Figs. 1 and 2 
we present the functions F^[z^ ^r) and F^^z^^r). 

The symmetry under the transition r 1/r is clearly visible. 
Also the conditions that ^ if either z ^ cxd or r ^ or 
r ^ 00 which follows from Eq. ( |154|) is evidently satisfied. For 



fixed z the functions oscillate with a frequency which grows with 
z. Since the amplitude decays rapidly, these oscillations are only 
visible in the log-plots. The correlations always decay like power 
laws, never exponentially. 
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Fig. 18. The function F^(z,r), linear (left) and |i<5(2;,r)| logarithmic are shown. 
Because of their small amplitude, the oscillations are invisible in the linear plot. To 
show the symmetry r — > 1/r, the r-axis is chosen logarithmic in both plots. 




The equal time correlation functions, to ^5(^^,1) are 

plotted in Fig. ^ left. On the right hand side we show Aa{r). All 
the functions, except F5 which is constrained by Eq. ( |148| ) pass 
through (For Fi the passage through is not visible on the plot 
since it occurs only at z ~ 30). The asymptotic behavior of the 
functions can be obtained analytically. The same is true for the 
functions Ai to A^. As argued above, all the functions Ai except 
A^ are symmetric under r 1/r. 
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Fig. 20. The functions |-Fj(z, 1)| are shown left. The zeros are visible as spikes in 
the log-plot (Further below, at z ~ 30, also Fi passes through zero.). Right, the 
functions |Aj(r)| = |i<i(0,r)| are shown. As discussed in the text, all of them except 
are symmetrical under r 1/r. 

6.1.1 Scalar, vector and tensor decomposition and CMB anisotropics 

The energy momentum tensor of seeds can be split into scalar, 
vector and tensor perturbations, since the time evolution of each 
of these components is independent. Furthermore, due to statisti- 
cal isotropy, the scalar, vector and tensor modes are uncorrelated. 



We use the decomposition of 6^-, given in Chapter 0, Eqs. 

The functions Fi to F5 determine the correlations: To work 
out the correlation functions we use also Eqs. (^,^ and ( |5T| , 



as well as ( |5^) . Using these identities and the ansatz (|T£7|), one 
easily verifies 



{U{r)<{r')) = {Mr)wKr')) = (/.(r)Tj.(rO) = (155) 
{fArH{r')) = {w,{r)r,,{r'))=0 (156) 

{fpir)f;ir')) = ^[2F5(.^ r) + 3F4(.^ r) 

+z\F3{z\ r)/r + F^iz^ l/r)r) + 

^z'F2{z\r) + ^z'F,{z\r)] (157) 
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{Mr)f:{r')) 



[3F5 + 4z^F2 + z^Fi] 



(158) 



-V^[^z^Fi{z\ r) + \f2{z\ r) + rFi{z\ l/r)] 
{w,{t)w]{t')) = 

^^[F, + z'F2]{k%-hk,) 

J — -F^\5ii5jm + 6im6ji — Sij6im + k [dijkikm + 

^im^ikj ^iikjkui S-iuikikj 6jik{kfji ^jmkiki^ ~\~ 
k kikjkikjji 



(159) 



(160) 



(161) 



It is interesting to note that although Cijim is analytic, the cor- 
relation functions of the scalar, vector and tensor components, in 
general, are not. The reason for that is that the projection op- 
erators onto these components are not analytic This is im- 
portant. It implies, e.g., that the anisotropic stresses in general 
have a white noise and not a k^ spectrum as erroneously con- 
cluded in [|iU5| , [73| . The scalar anisotropic stress potential thus 
diverges on large scales, (I/ttP) oc l/(r/c'^) for /cr <C 1. A re- 
sult which is also obtained in the large- limit and in numerical 
simulations of 0{N) models. The power spectrum of the scalar 
anisotropic stress potential fj^ is analytic if vector and tensor 
perturbations are absent, F^ = F2 = 0. In the generic situation, 
F^iz = 0, r = 1) = ^5(1) 7^ 0. ^ ^"^ 



Another situation where /jr has a white noise spectrum is the 
case of perfectly coherent seeds [)56j , in other words if 



Cijim{k., r, t') — Aij{k, r)A^^{k, t'] 



(162) 



^^^Mjn though the potential and thus also the Bardeen potential ^' (see Eq. 
( |166| )) diverge for fcr — > 0, the physically relevant (measurable) quantities like Q^y 
and R^y stay perfectly finite. 
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The fact that as weU {fpf*) as {fpfl) are white noise imphes that 
also {fnfn) = M{fpfn)\^/{fpfp) ^lust behave hke white noise 
and thus F5 oc and F2 oc 2^ on large scales. This can also 
be obtained by using the analytic properties of the correlators 
(Oo^O*^) and energy momentum conservation . 



Generically, we expect the following relation between scalar, vec- 
tor and tensor perturbations of the gravitational field on super- 
horizon scales, X = kr <^ 1: (The equations for the scalar, vector 
and tensor gravitational potentials in terms of /., w and are 
given in Chapter 0.) 



(163) 

(|S,|^>~— A5(l) (164) 
(|HyP>~4fV>l5(l) . (165) 



If the large scale CMB anisotropics were solely induced by super 
horizon perturbations, this could be translated into a ratio be- 
tween the scalar, vector and tensor contributions to the C^'s on 
large scales, i ~ 50. However, since the main contribution to the 
CMB anisotropics is induced at horizon crossing, x = 1 (see be- 
low) the above relations cannot be translated directly and we can 
just learn that one expects, in general, contributions of the same 
order of magnitude from scalar, vector and tensor perturbations. 

We want to discuss in some detail the CMB anisotropics induced 
from scalar perturbations. From 

$ + ^ = -2eA , (166) 



we see that even if $ has a white noise spectrum due to 'compen- 
sation' this generically leads to a k^"^ spectrum for ^ and 
for the combination $ — ^ which enters in Eq. ( |167|) below. 



This finding is in contradiction with [|i0b|,[73|j , which predict a 
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white noise spectrum for ^, but it is not in conflict with the Har- 
rison Zel'dovich spectrum of CMB fluctuations which has been 
obtained numerically in [|iiY| , |5Tl ,^] . This can be seen by the follow- 
ing simple argument: Since topological defects decay inside the 
horizon, the Bardeen potentials on sub-horizon scales are dom- 
inated by the contribution from dark matter and thus roughly 
constant. The integrated Sachs Wolfe term then contributes only 
up to horizon scales. Therefore, using the fact that for defect 
models Dg and V are smaller than the Bardeen potentials on 
super-horizon scales, we obtain 

1 

+ / - *')(^) x)ji{xQ - x)dx , (167) 

where x = kr and prime stands for derivative w.r.t. x. The lower 
boundary of the integrated term roughly cancels the ordinary 
Sachs Wolfe contribution and the upper boundary leads, to 



k'{\{AT/TUk)\%sw 
6^(3^5(1) + 4^2(1) + F3{l)]jf (xo), (168) 

a Harrison-Zel'dovich spectrum. The main ingredients for this 
result are the decay of the sources inside the horizon as well as 
scaling, the rest follows for purely dimensional reasons. 

The parameter space of generic causal scalar seed models pro- 
vided by the five functions Fi to F5 (of two variables) is still 
enormous and is rather impossible to investigate. For a realistic 
model, the parameter space is even larger due to the radiation- 
matter transition which breaks scale invariance: the seed func- 
tions can be different in the radiation and in the matter era. For 
global 0{N) defects this difference turns out not to be very im- 
portant (less than about 20% [g7[]) it may, however, go to factors 



of two and more for cosmic strings ||i5^| . 
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6.2 Mimicking inflation 



Neil Turok has constructed a model with scaling causal seeds 
which perfectly reproduces the CMB anisotropy spectrum of in- 
flationary models [p7H|. Other synthesized causal seed models 



with various heights of the acoustic peaks are discussed in [^191 , ^16 



Spergel & Zaldarriaga argued that causal seeds can nevertheless 
be distinguished from inflationary models by the induced polar- 
ization [|T35|. In their argument they however use that the corre- 



lator (/tt/tt) be white noise. As we have seen, this is only correct 
for purely scalar or perfectly coherent seeds. Therefore, allowing 
for vector and tensor contributions, as well as for decoherence, 
one can in principle circumvent the Spergel & Zaldarriaga argu- 
ment. However, the fact that seeds notoriously lead to too low 
amplitude acoustic peaks, limits very strongly the allowed vec- 
tor and tensor contributions which enhance the CMB anisotropy 
spectrum in the Sachs- Wolfe, but not in the acoustic peak re- 
gion. Furthermore, the smearing out of the acoustic peaks in- 
duced by decoherence is not observed in the data, limiting the 
allowed amount of decoherence considerably. 

We restrict the following discussion to 'perfectly coherent' models 
with purely scalar perturbations. We parameterize the seeds by 



r)^:(fc, rO) = ^^,Pi{z, r) (169) 

($,(fc, r)$:(fc, t')) = y^^2(^, r) (170) 

(*,(fc, r)$:(fc, t')) = -^j^^Psiz. t) . (171) 
Perfect coherence then implies 



Pi{z, r) = ±V^i(v^, l)P2(y^, 1) (172) 
Due to the absence of vector and tensor perturbations, the sum 
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$ + oc /tt is suppressed by a factor on large scales, z <C 
1 [153]. In a first attempt we simply set = — which implies 



Pl = P2 = -Ps = P. 

We discuss two famihes of models (see 
Family I 

To enhance the acoustic peak, we use seeds which are larger in 
the radiation era than in the matter era. 



where here the subscripts r and ^ indicate the radiation and mat- 
ter era respectively. The parameters t and b are varied to obtain 
the best fit and the amplitude e is determined by the overall nor- 
malization. 
Family II 

The second family of models is inspired by Ref. [ |i44|j , which stud- 
ies spherical exploding shells with p + 3j9 oc 6{r — At). In terms 
of the source functions defined in Chapter |^, this model is char- 
acterized by 



1 sm{AkT) 

Q;^1/2 J^]^^ 

E(t 



^2^3/2 C2 



cos{Ckr] 



sin(C/cr) 
Ckr 



with a = (d/a)r and E = {4 — 2/a)/{3 — 12a). The functions 
fp and are then determined by energy momentum conserva- 
tion, Eqs. ( ^ , |5I3| ) . The function E is chosen such that the power 
spectrum of fj,- is white noise on super horizon scales, a condition 
which is required for purely scalar causal seeds as we have seen 
in the previous section. This leads to the Bardeen potentials 
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* = -$-2e/, . 



(175) 
(176) 



Here the seed functions are actually not given as random variables 
but as square-roots of power spectra, and one has always to keep 
in mind that we assume perfect coherence. Of course one can also 
regard Eqs. (|i 7i3|,|i V6|) as mere definitions with 



Pi(z,l)=rA:%Ve', (177) 
P2{z,l)=Tk^^^/e\ (178) 

P^{z, l) = Tk^'^^/^ = -v/Pi(2,l)P2(^,l) . (179) 



With a somewhat lengthy calculation one can verify that E is 
chosen such that f-j^ oc const, for z <C 1 and the functions Pi{z,l) 
are analytic in z"^ = (kr)'^. This family of models is described by 
the parameters A and C, which have to satisfy < ^4, C < 1 for 
causality. Also here one can choose different amplitudes for the 
source functions in the radiation and matter era by introduction 
of the additional parameter t 1. 



Seeds generically produce isocurvature perturbations. For a flat 
universe, this implies a position of the first peak at ^ ~ 350, which 
is definitely incompatible with the recent CMB observations (see 
also [ |i2U|] , [^). However, the tight constraints on the flatness 
of the universe obtained from CMB data analysis are based on 
the assumption of adiabatic primordial fluctuations. Using this 
loophole, it is possible to construct closed A-dominated isocur- 
vature models which have the first acoustic peak in the observed 
position. 



For a given seed-model, the position of the first acoustic peak 
is determined primarily by the angle subtended by the acoustic 
horizon Xac at decoupling time, r^ec- The angle under which a 
given comoving scale A at conformal time Tdec is seen on the sky 
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is given by ^(A) = A/x(ro - r^ec), where 



x{y) = 



sm{y) if K > 
smh{y) if K < 
y if K = 



is the comoving angular diameter distance (K denotes the cur- 
vature of 3-space). As the harmonic number i is inversely pro- 
portional to the angle 6, this yields £pe&k — R'^ltlk where R = 
O^f /6ac and ^pl*]^ is the peak position in a flat model with the 
same value of ^mh^- In terms of cosmological parameters one 
finds (see Ref. 
of [|Tg), 



or a factor 1/2 is missing in the formula 



R 



1 

2^ 



K\ 



An interesting point is that for the quantity R depends 

very sensitively on Qa- Thus, we can have important shifts in the 
power spectrum, R ~ 0.6 say, with relatively small deviations 
from flatness {^^ = 0.3, Qa = 0.9, = -0.2). In Ref. \^ the 



authors have shown that the simple prescription £ Ri repro- 
duces the CMB power spectra for curved universes within a few 
percent. On lines of constant R, CMB power spectra are nearly 
degenerate. This simple prescription can be used to rescale the 
flat spectrum. Furthermore, one has to ascertain that the value 
of h^Qrn used in the spectrum calculation agrees roughly with the 
value preferred by the best fit value of R and the super-novae con- 
straint [|i25|] , which can be cast in the form ~ O.TS^a — 0.25. 
h'^flrn determines the time of equal matter and radiation and 
thus influences the early integrated Sachs- Wolfe effect, which con- 
tributes to the spectrum right in the region of the first peak. We 
therefore get a better approximation if we use the correct value 
for Qjnh'^. 
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The results from an analysis of family I are shown in Figs. |2T] 
and |22| (long dashed hnes). The best fit model is shown. Clearly, 
even though we can change the cosmological parameters to fit the 
first acoustic peak, the second and third peaks are also shifted by 
this procedure and can no longer be fitted. The cosmological pa- 
rameter R is now fixed by CMB anisotropy measurements due to 
the measured inter-peak distance. This situation was still slightly 
different with the 'old' Boomerang and Maxima data [|35| , |68|] for 



which this model could still provide a reasonable fit in ad- 



dition, the shifted first peak is too narrow to fit the data well. 

The 'best fit' model corresponds to the parameters t = 2.2, b = 
1/9, = 0.35 and R = 0.53. It has a value of = 96, which, 
for 13 points and 4 parameters (t, 6, R and the normalization), 
excludes it at more than 99.8% c.l. if Gaussian statistics are as- 
sumed. 



250 I I ' I I I I I I I I I I ' I I I I ' I 




Q I I I 1 I I I 1 I I 1 I I I I I I I I ± 
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I 

(T) 

Fig. 21. The CMB temperature anisotropy spectrum £{£ + l)Cf, for our best fit 
model of family I (long dashed) and family II (solid) is compared with the new B98 
data. Family I model is closed, O ~ 1.2 it has problems fitting anything. The family 
II model is flat and is in good agreement with the data (x^ = 4/9), we fitted only 
i < 700. A standard inflationary spectrum with the same cosmological parameters 
as the family II model {h = 0.65, h'^Qb = 0.019, ^cdm = 0.35, Qa = 1 — ^m) is also 
indicated (dotted). 
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200 400 600 800 1000 
I 

( F) 

Fig. 22. The CMB polarization spectrum s for our best fit model of family I 
(long dashed) and family II (solid) is compared with a standard inflationary spec- 
trum with the same parameters as above (dotted). The family I model predicts a 
larger polarization signal in the band 50 < £ < 500. On the contrary, the lack of 
intermediate scale polarization at ^ < 200 in the family II model is clearly visible. 

A much better fit can be achieved by the models of family II. 
The best fit model of this family has = 4.0 for 13 points 
and 5 parameters (^4, C, R and the normalization) and is in 
good agreement with the data (see Fig. sohd hne). The model 
shown corresponds to the best fit parameters A = 0.9, C = 
0.8, t = 0.9, and R = 1. In this model which is flat and causal, 
the first peak in the polarization spectrum is suppressed, as has 
been noted in Ref. [|i35|] (see Fig. [22], solid line). 

We can therefore conclude that generic causal scahng seed models 
for structure formation can reproduce the recent CMB anisotropy 
data [|351 , |55|] . To achieve this agreement, we have suppressed vec- 



tor and tensor perturbations and have assumed perfectly coherent 
fluctuations. We believe that it is quite improbable that topolog- 
ical defects from a GUT phase transition have such a behavior. 
Nevertheless, there might be some other scale-invariant causal 
physical mechanism (e.g. some spherically symmetric 'neutrino 
explosions', see Ref. [[T^j ) leading to seeds of this or similar type. 
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Clearly, we only have a satisfactory model of structure formation 
if also the physical origin of the 'seeds' is clarified. 



6.3 Mixed models 



In this section we want to discuss the possibility of combining 
topological defects with primordial, inflationary perturbations. 
There are at least two important reasons why such an analysis 
is interesting. First of all, as we saw in the previous chapters, 
the predictions based on topological defect models for the CMB 
power spectra are dramatically different from what is expected in 
the standard inflationary model. This means that, in principle, 
the two contributions can be disentangled. Secondly, this analysis 
contributes in a phenomenological but physically motivated way 
to the discussion, about the model dependence of the values of 
cosmological parameters derived from accurate CMB data. Infla- 
tionary models leading to a mixture of inflationary perturbations 
and defects which both contribute to the CMB anisotropics have 
been developed e.g. in Refs. [j7^,^,^ and others. A similar 



analysis, but with cosmic string and with an older CMB dataset, 
can be found in []27[] . A mixture of global defects and inflationary 



perturbations has also been considered in [jT^ 



We consider an initial fluctuation spectrum with scalar initial per- 
turbations given by their spectral index ng, in addition to topo- 
logical defects which formed at some symmetry breaking tem- 
perature Tc determining the amplitude of the seed perturbations. 
For definiteness we take the defects to be cosmic texture. It is 
reasonable to assume that the inflationary fluctuations and those 
induced by topological defects are not correlated and therefore 
the resulting perturbation spectra can be added in quadrature. 
We set 



l + r) ' (1 + r) 



^Cfo^^^ (180) 
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for the CMB anisotropy, where A is the pre factor in units of 
(jCOBE^ C'/o = = 1 and the parameter r gives the relative 
amphtude of textures/inflation at ^ = 10. We let vary the cosmo- 
logical parameters as follows: = 0.015 to 1.0; Qb = 0.015 to 
0.2; = 0.0 to 1.0 and h = 0.25 to 0.95. We restrict the analysis 
to flat universes. We vary the spectral index of the primordial in- 
flationary perturbations within the range Ug = 0.50 to 1.50. The 
theoretical inflationary models are computed using the publicly 
available CMBFAST program and are compared with the recent 
BOOMERanG-98 [TTOj, DASI ^ and MAXIMA-1 results. 
The power spectra from these experiments were estimated in 19, 
12 and 10 bins respectively, spanning the range 25 < £ < 1000. 
For the DASI and MAXIMA-I experiments, we use the pub- 
lic available correlation matrices and window functions. For the 
BOOMERanG experiment, we assign a flat window function in 
each bin £{£+ l)Q/27r = Cb, we approximate the signal Cb to be 
a Gaussian variable, and we consider ~ 10% correlations between 
neighboring bins. The likelihood for a given cosmological model 
is then defined by -21nL = (Cf - Cf)MBB'{Ci, - Cf,) where 
Mbb' is the Gaussian curvature of the likelihood matrix at the 
peak. We consider 10%, 4% and 5% Gaussian distributed cah- 
bration errors for the BOOMERanG-98, DASI and MAXIMA-1 
experiments respectively. We also include the COBE data using 
Lloyd Knox's RADPack packages [ [Ibb . 



In Fig. |23| we plot the obtained likelihood distribution for the 
parameter i?, after marginalization over the remaining 'nuisance' 
parameters and as function of different external priors. The present 
CMB data plus an external Gaussian prior on the value of the 
Hubble constant h = 0.65 ± 0.2 gives r < 2.6 at 95% c.l.. Hence, 
we cannot exclude a substantial admixture of topological de- 
fects from CMB data alone. Adding strong constraints on the 
baryon density parameter from big bang nucleosynthesis Q^bh? = 
0. 020±0. 004 gives r < 2.1. Finally, combining the CMB data with 
LSS observation by including a constraint crgO^^ = 0.50 ± 0.05 
gives r < 1.6, always at 95% c.l.. As we can see, even assuming 
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quite restrictive priors, the present CMB data allow a contribu- 
tion from textures as big as r ~ 1. Furthermore, models with 
non-zero r are in slightly better agreement with the data, but 
there is no significant evidence for non-vanishing r. 

The main result of this analysis is that current microwave back- 
ground data do not exclude a dominant contribution from tex- 
tures on large scales, and marginally favor a significant fraction. 

It is interesting to study how textures can affect the constraints 
on the remaining cosmological parameters. Marginalizing over 
Qm, Qa, h and Qb- Assuming the external priors h = 0.65 ± 0.2 
and Qfo/i^ = 0.020 ± 0.002, we obtain the likelihood contours 
shown in Fig. ^ (left panel) in the r- Ug parameter space. There 
is a degeneracy between the amplitude of the texture component 
and the spectral tilt. A blue tilt (n^ > 1) is found to be compatible 
with a larger textures contribution. 

If we instead marginalize over ^^a, ^5, without assuming ex- 
ternal priors from big bang nucleosynthesis, we find the contours 
shown in Fig. |2?| (right panel) in the r - Qbh'^ plane. Without 
the BBN constraint, the contribution of textures modes can be 
even larger, and there is clearly a degeneracy along the r — Q^/i^ 
direction. 

The best fit spectrum with = 0.175, Qa = 0.825, h = 
0.75, Qb = 0.045,77,5 = 1.02, r = 0.9 together with a 'concor- 
dance' fit, = 0.3, Qa = 0.7, h = 0.65, Q5 = 0.04, n, = 1, r = 
are shown in Fig. 



100 





Weak Priors 




-^t;; — ^ Strong BBN 




~\ \. 0.45 <a,Q/ '< 0.55. 











I ' 1 ' 1 ' 1 ' I 

0,0 0,5 1,0 1,5 2,0 



Textures contribution r 

Fig. 23. The likelihood probability distribution functions for r in the case of different 
external priors. 
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Fig. 24. Left: the likelihood contours for r and Ug are shown after marginalization 
over Om,) 0^ and Oa with weak priors on h and $7;, (see text). Right: the likelihood 
contours for r and Uh = Qbh'^ are shown after marginalization over Ug, 0,^, h and 
i^A. No priors are applied 
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Fig. 25. The best fit mixed spectrum is shown. The corresponding parameters are 
= 0.175, Qa = 0.825, h = 0.75, Qb = 0.045, = 1.02, r = 0.9. Its inflationary 
(long dashes) and texture (short dashes) components as well as the spectrum from 
a concordance model (dotted, see text) are also indicated. 
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7 Conclusion 



In this review we have discussed the role of topological defects 
for cosmic structure formation. Even though the basic idea seems 
natural and intriguing, we have seen that structure formation 
with global 0(n) defects and (local) cosmic strings is ruled out 
already by the CMB anisotropy data alone. Due to the isocur- 
vature nature of the perturbations and due to the importance of 
vector and tensor perturbations, defects do not induce the acous- 
tic peaks visible in present data. Furthermore, non-linearities in 
the defect evolution tend to decohere, smear out, the peak struc- 
ture. Nevertheless, generic causal seeds, which are dominated by 
scalar perturbations and have nearly no decoherence, may mimic 
the CMB anisotropy spectrum of inflationary models. Such mod- 
els, however, do certainly not represent topological defects for 
which non-linearity and therefore decoherence is one of the ba- 
sic ingredients. Furthermore, such models could be distinguished 
by their polarization spectrum. Models where both, inflation and 
topological defects play a role for structure formation, such that 
the large scale CMB anisotropics from both components are of 
comparable amplitude, cannot be excluded with present data. 
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APPENDIX 



A Definitions of all gauge-invariant perturbation variables 

In this Appendix we give precise definitions of all the gauge- 
invariant perturbation variables used in the text. These defini- 
tions, their geometrical interpretation and a short derivation of 
the perturbation equations can be found in Refs. [j?2],^. We re- 



strict the analysis to the spatially fiat case, K = 0. The general 
case can be found in the references above. We define the per- 
turbed metric by 

g = g + a^h, (A.l) 

where g denotes the standard Friedman background, a is the scale 
factor and h denotes the metric perturbation. 



A . 1 Scalar perturbations 

Scalar perturbations of the metric are of the form 



hiS) ^ _2A{drf + 2iBkjdTdx^ + 2{Hl + -HT)6,jdx'dx^ 

3 

-2k-^HTk,kjdx'dx^ . (A.2) 

Computing the perturbation of the Ricci curvature scalar and the 
shear of the equal time slices, we obtain 



6R = 4a-^k^n with n = HL + \ht (A.3) 
S = aa{^ - -6i)dx' ® dj, with a = -Ht - B . (A.4) 
The Bardeen potentials are the combinations 
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^ = A - k~\{a/a)a - a] 



(A.5) 
(A.6) 



They are invariant under infinitesimal coordinate transformations 
(gauge transformations). In longitudinal gauge defined hy B = 
Ht = 0, the variable $ represents the perturbation of the spatial 
metric while is the perturbation of the dt^ term, the lapse 
function. 

To define perturbations of the most general energy momentum 
tensor, we introduce the energy density p and the energy flux u 
as the time-like eigenvalue and normalized eigenvector of T^, 

TV = -pu" , u^ = -l. 



We then define the perturbations in the energy density and en- 
ergy flux field by 

P = pa + 6) , (A.7) 
u = u^dt + u% ; (A.8) 



is fixed by the normalization condition, = a~^{l — A) and 
p is the homogeneous background energy density.. In the 3-space 
orthogonal to u we define the stress tensor by 

T^" = Pf'^P'^^T''^ , (A.9) 



where P = u<S)u + g is the projection onto the sub-space of TA4 
normal to u. It is 

^0 _ _ i _ Q 



The perturbations of pressure and anisotropic stresses can be 
parameterized by 

r/ = p[{l + 7rL)S-' + tt/'] , with 4 = 0, (A.IO) 
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where p is the background pressure. For scalar perturbations we 
set 

k 

and 

ii^f = {-k^''k'kj + \5'^)Il . 

Studying the behavior of these variables under gauge transforma- 
tions, one finds that the anisotropic stress potential 11 is gauge 
invariant. Another gauge invariant combination from the matter 
variables alone is 

T = TTL--5. (A.ll) 

w 

r is proportional to the divergence of the entropy flux and van- 
ishes for adiabatic perturbations []53]. A gauge invariant velocity 
variable is the shear of the velocity field, 

= {k-%kj - ^%)aV , with V = v- k'^HT. {A.12) 
o 

There are several different useful choices of gauge invariant den- 
sity perturbation variables, 

L», = (5 + 3(1 + w){d/a)k'^a (A.13) 
Dg = 6 + 3(1 + w)n = Ds + 3{l + w)<P (A.14) 
D = D, + 3(1 + w){d/a)k-^V . (A.15) 

In this work we mainly use Dg. Here w = p/ p denotes the en- 
thalpy. Clearly, these matter variables can be defined for each 
matter component separately. For ideal fluids like CDM or the 
baryon-photon fluid long before decoupling, anisotropic stresses 
vanish and ttl = {cl/w)5, where Cg is the adiabatic sound speed. 

Also scalar perturbations of the photon brightness, t^^^ are not 
gauge invariant. To define the brightness we first consider the one 
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particle photon distribution function 

f{x, p) = f{r, p) + F(r, X, p) , 

where / is just the Bose-Einstein distribution. The background 
Liouville equation requires that / be a function of the redshift 
corrected momentum v = ap only. The brightness perturbation is 
then the momentum integral of F, which depends on the photon 
direction, position and time 




(A.16) 



This can be decomposed in scalar vector and tensor contributions, 



It has been shown [|52|1 that the combination 



^is) ^ ^{S) + + 4ik'^n^kja (A. 17) 

is gauge invariant. This is the variable which we use here. In other 
work[|7^ the gauge invariant variable 6 = + $ has been used. 
Since $ is independent of the photon direction n this difference 
in the definition shows up only in the monopole, Cq. But clearly, 
as can be seen from Eq. (|A.1Y| ), also the dipole of Ci, is 
gauge dependent. 

The brightness perturbation of the neutrinos is defined the same 
way and is not repeated here. 

A. 2 Vector perturbations 

Vector perturbations of the metric are of the form 

/i(^) = 2Bjdx^dt + ik-\kiHj + kjHi)dx^dx^ , (A. 18) 
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where B and H are transverse vector fields. The simplest gauge 
invariant variable describing the two vectorial degrees of freedom 
of metric perturbations is XI, 

^, = k-'H,-Bj. (A. 19) 

Vector anisotropic stresses are gauge invariant. They are of the 
form 

Trip = ik-\kjUi + kUj) . (A.20) 

The vector degrees of freedom of the velocity field are cast in the 
vorticity 

ui-j — Uj-i = ia{kjUJi — kiUj) with uj = Vj — Bj . (A. 21) 

Vector perturbations of the photon brightness are gauge-invariant. 
To maintain the notation consistent we denote them by Ai '-^^ = 
.(^). 

A. 3 Tensor perturbations 

We define tensor perturbations of the metric by 

h^^^ = 2Hijdx'dx^ , (A.22) 

where Hij is a transverse traceless tensor field. 

The only tensor perturbations of the energy momentum tensor 
are anisotropic stresses, 

4"^^ = n,,- . (A.23) 

Again for notational consistency tensor perturbations of the pho- 
ton brightness are denoted M.^^^ = 
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Clearly, all tensor perturbations are gauge-invariant (there are no 
tensor type gauge transformations). 



B Boltzmann equation and polarization 



The relativistic Boltzmann equation for the photon distribution 
function is of the form 



:B.r 



where f{t,x,p) is the one-particle distribution function on the 
mass bundle Pm = {{p,x) G TM\g{x){p,p) = —m^} and C[f] is 
the collision integral which describes interactions. The left hand 
side of (|B.i| ) requires the particles to move along geodesies in the 
absence of collisions. (For a thorough treatment of the kinetic 



approach in general relativity see e.g. [|i39 
therein.) 



or 



i4Y| and references 



B.l The collisionless case, momentum integrals 



Let US first consider the situation where collisions are negligi- 
ble, C[f] = 0. The unperturbed Boltzmann equation then im- 
plies that / be a function oi v = ap only. Setting / = f{v) + 
F{T,x,v,n), where n denotes the momentum directions, leads 
to the perturbation equation 



drF - n'd^F = WA 



[B 



H,j) + H 



L 



;b.2) 



To derive ( [b.2| ) we have used p^ = 0. The Liouville equation for 
particles with non- vanishing mass can be found in Ref. [j?^ . 
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The ansatz 



T[x, n) I \T{x, n] 



with T{x,n) = T{t) + AT(a:, n), and where g^^'' is the spatial 
part of the metric, leads to 

Note that in terms of the brightness perturbation, i defined in 
Appendix^, 

4:7T Ar 1 

I = —7 / Fv^dv we have = -i . (B.5) 

pa^ ^ r 4 

Comparing this with equation ( [b.^|) , we obtain 



- n^i, - - H,,) n'n^ + Uj\ ■ (B.6) 

The fact that gravitational perturbations of Liouville's equation 
can be cast entirely into temperature perturbations of the original 
distribution is not astonishing. This is just an expression of grav- 
ity being "achromatic", i.e. independent of the photon energy. 

We now decompose into scalar, vector and tensor compo- 

nents. Since components with spin higher than 2 are not sourced 
by the right hand side of equation (|B.(3| ) and since they are sup- 



pressed at early times, when collisions are important, we neglect 
them. 



For the scalar contribution to AT/T we obtain from ( [b.5| ) (in 
Fourier space) 
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1 



dr [-y) + ^^ ^^ ) 



;b.7) 



This equation is not manifestly gauge-invariant. However, defin- 
ing 

^M^'^=[^) +H, + ^eH + in%{H-B), (B.8) 
it reduces to 

drM^^^ + fikM^^^ = 4i/i/c($-^), (B.9) 

where $ and ^ are the Bardeen potentials and we have set fi = 
n^kj/k. Since the right hand side of ( |B.9|) is gauge invariant, the 
left hand side must be so as well and we conclude that Ai''^^ is 
a gauge-invariant variable (a direct proof of this, analyzing the 
gauge transformation properties of the distribution function, can 
be found in \W^)- 



lji4(^) coincides with the scalar temperature fluctuations up a to 
a gauge dependent monopole and dipole contribution. 

The vector and tensor contributions to AT/T are gauge invariant 
by themselves and we denote them by ^A4^^^ and ^A4^^\ In the 
absence of collisions, they satisfy the equations 



M'^^^ + iiikM^^^ = -iiiin'^'k^^^^ (B.IO) 
M^'^^ + ilikM^^^ = in^n'^H^e. (B.ll) 

The components of the energy momentum tensor are obtained 
by integrating the second moments of the distribution function 
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over the mass shell, 



(B.12) 



with the invariant measure ijl{x^p) = ^^^^^^^Sp. One finds 



47r 
3z 



(0) 

■> 



IGtt 

-9 



yin) 

nil) 

Stt 

r(7) = o , 
1 



n(^) = 



1^ 1 iinM^^UVt 
—6i 



5 
V6 



+ CT]^ 3 , (72^1 + (72,3 ? J 



ttST^ Tri^) 




(B.13) 

(B.14) 

(B.15) 
(B.16) 
(B.17) 



(B.18) 



(B.19) 



(7ag) + 104? + ^cjI 



with .r=- 

and .r = ^(M^nMr^ 



(T) 
i 
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Here, we have also expressed the result in terms of the multipole 
moments defined by the expansions (|B.t3(J| ) and (|b.52| ) for a wave 
vector k pointing in 2^-direction.. 

The expressions for the neutrino fiuid perturbations in terms of 
M are identical. 



B.2 Collisions, polarization 



Let us now turn to the collision term. Just before the process 
of recombination during which the fiuid description of radiation 
breaks down, the temperature is ~ 0.4 eV. The electrons and 
nuclei are non-relativistic and the dominant collision process is 
non-relativistic Thomson scattering. 

Thomson scattering depends also on the polarization of the in- 
coming radiation field. It is therefore necessary to treat polar- 
ization. This is usually done by introducing the Stokes parame- 
ters [|781jg3|,PI7|,^]: 



For a harmonic electro-magnetic wave with associated electric 
field 

E(^r, t) = (eiEi + €2^2) e*^^ ' ^"^'^^ , (B.20) 



where n, ci and 62 form an orthonormal basis and the complex 



field amplitudes are parameterized as Ej = aje'^^^j the Stokes 



parameters are given by 

I = al + al (B.21) 

Q = al-al (B.22) 

[/ = 2aia2Cos(^2-^i) (B.23) 

V = 2aia2sin((52-^i). (B.24) 



114 



/ is the intensity of the wave (whose perturbation A4 has inter- 
ested us so far), while Q is a measure of the strength of hnear 
polarization (the ratio of the principal axis of the polarization 
ellipse). U and V give phase information (the orientation of the 
ellipse). V also gives the amplitude of circular polarization. For 
non-relativistic Thomson scattering V is completely decoupled 
and (since it vanishes at early times) is therefore never gener- 
ated. 

As Q and U vanish in the background, perturbations cannot cou- 
ple to them (since such terms are 2nd order), and the equation 
corresponding to (|b.Y| ) for Q and U become (neglecting scatter- 
ing) 

driQ; U) + in^h{Q] U) = 0. (B.25) 



The differential cross section of Thomson scattering for a pho- 
ton with incident polarization e^j) scattering into the outgoing 
polarization e(s) = e' is [j7B] 



da 



Stt 



;b.26) 



n 




7 



Fig. B.l. Definition of the angles and vectors for Thomson scattering in the (n, £2) 
plane. 

It is often convenient to introduce the two 'partial' intensities 
Ii = al = {I + Q)/2 and /2 = a2 = {I — Q)/2. A wave scattered in 
the (n, €2) plane (see figure |B.i|) by an angle 6 has the intensities 



3aT 



I 



(0 



Stt 

3crr Jt) 2 

R ' cos 

Stt ^ 



(B.27) 
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or, expressed in terms of the Stokes parameters, 



( /(•) ^ 


3aT 




IGtt 

V 



^ 1 + cos^ 6 sin^ 6 ^ 
sin^ 9 1 + cos^ 6 



B.28) 



A rotation in the (ei,e2) plane doesn't change the intensity of 
the wave, but it changes Q and IJ to 



g' = Qcos(20) + t/sin(2(/)) 
[/' = -t/sin(20) + Qcos(20) 



(B.29) 
(B.30) 



z 


\w'-(p 








\a \ 






^6*' / \ 
/n \ 

e /\ 1 















Fig. B.2. Definition of the angles and vectors for Thomson scattering in the general 



case. The polarization vectors are oriented like in figure B.l. 



To determine the cross section that a given 'initial' wave 
{I^^\Q''^\U^^^) propagating in direction n be scattered into a 
wave Q("), with direction n\ we need to go through the 
following steps (we will use the plane (?/, z) as reference plane, 
see figure ( |B.2| ) for definitions of the angles and vectors): 

(1) Rotate around n such that the plane {n,n') turns into the 
plane (nz). One needs to apply the rotation ( [b.29| , [BT^) for 
(j) = a to the Stokes parameters. 

(2) Rotate the new plane (n, n') around z into the reference 
plane (2/, z). This operation does not influence the incoming 
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Stokes parameters.. 

(3) Now we are in the known case of ( [b.l^Y|) and Hence 
we can apply the scattering matrix. 

(4) We then rotate the scattering plane back around z into the 
old {z,n') plane. This does not change the scattered Stokes 
parameters. 

(5) Finally we rotate around n' by the angle a' to reach the orig- 
inal state. To do this, we have to apply the rotation matrix 
(|BT2^,|BT33) again, but now for = a'. 



Following the steps outhned above, we recover the scattering ma- 
trix in the basis (/i, I2, U) given in equations ( [b.31^|) - ( [b.35|) (see 
also 123]). V is completely decoupled from the other parameters 
and follows an evolution which is independent of the rest. Hence 
by starting with V{t <^ tdec) = it will stay zero and can be 
neglected. The scattering matrix P, which determines the (non 
vanishing) scattered Stokes parameters from the initial ones. 



/ As) \ 
^1 




1 \ 
^1 


7-(^) 


= 

4:71 











is then given by 



P 



;b.32) 



where 



p(o) = 



3 
4 



2(1 



/i2)(l-/i 





/2> 



1 
Oy 



;b.33) 
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p(l) = - 



4;U/i'cos(0' — 0) 2;usin(0' — 0) 



-Vsin((/)' - 0) 2cos((/)' - 0) 



p(2) = 



/ At^/i'2cos[2(</)' -0)] -/i2cos[2((/)' -</))] /i2^'sin[2(</,' - 0)] 
cos[2{0' - (/))] cos[2((^' - fli)] -^i' sm[2((^' - fli)] 
~2AtAt'2 sm[2((^' - 0)] 2^fsm[2{0' - (/))] 2(1(1' cos[2{(j>' - ^)] 



(B.34) 



;b.35) 



As we are in an isotropic situation, we will perform all the calcu- 
lations in a special coordinate system with k \\ z and n, n' as in 
Fig. IB.4 Clearly the results are independent of this coordinate 
choice. 

The matrix R connecting (/i, /2, U) with (/, Q, U) is given by 



;b.36) 







'l/2(/ + Q)' 


1 

~ 2 


fl 


10^ 






(A 


h 




1/2(7 -Q) 


1 


-1 


Q 


= R 


Q 












2^ 









To calculate the collision term including polarization , we change 
into the {h.h) basis. For each of the two intensities A G {1,2} 
we then have the collision term given by 



C[f 



(A) 



df- 



(A) 



dr 



dr 



(B.37) 



where and f'^^ denote the distribution of photons in the 
polarization state A scattered into respectively out of the beam 
due to Thomson scattering. 

In the matter (baryon/electron) rest frame, which we indicate by 
a prime, we know that 



(A)/ 



-b,n) = ^/ P\{n,n')f^'^'{p',n')d^' , 



dt' 



47r 
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where denotes the electron number density and is the 2x2 
upper left corner of the normalized Thomson scattering matrix 
(|B.32|) . In the baryon rest frame which moves with four velocity 
w, the photon energy is given by 

p' = p^u^ . 



We denote the photon energy with respect to a tetrad adapted 
to the slicing of space-time into {r = constant} hyper-surfaces 
by p : 

p = p^n^" , with n = a-^[{l- A)dr + B'd,] , 



The lapse function and the shift vector of the slicing are given by 
a = a{l + A) and f3 = —B^di . In first order, 

Pq = ap{l + ^4) — apriiB^ , 



and to zeroth order pi = aprii. Furthermore, the baryon four 
velocity has the form = a~^{l — A) , = w^-u* up to first 
order. This yields 

p' = Pi,u^' = p{l + n^{v' - B')) . 



Using this identity and performing the integration over photon 
energies, we obtain 



Pr 



dr' 



47r 



1 + 4n,{v' - B')^ 
I L^^\n')Ps\n,n')dn' . 



1 



Photons which are scattered leave the beam with the probability 
given by the Thomson cross section (see e.g. 



TUT 



(A) 



arnj^'^'ip^n) , 



dt' 



119 



so that we finally have 



dfr df. 
dt' 



dt' 



+ L^^\n')Ps\n,n')dn'] 



(B.38) 



By setting C^-^) = C^^^+C^^) and C^^) = C(i)-C(2) we transform 
the collision integral back to the normal stokes parameters. The 
term for U has the same form as the one for Q, just with the 
corresponding matrix elements in the integral. The Boltzmann 
equation then finally becomes (setting 8 = {A4, Q,U) and for 
the flat case, k = 0): 



-M - AifiVb + 4n^LJb,i + / dn'p^e'^ 



Q + i/ikQ = aarUe -Q + J d^'P^S'^ 
if + iiJ,kU = aarUe [-U + J dfl'P^S'^ 



(B.39) 
(B.40) 
(B.41) 



where we have to use the scattering matrix transformed into the 
{M, Q, U) basis, 



P = Ps + Pv + P7 



(B.42) 



with 



( 



3 - - ^'2 ^ 3^ V2 (l-3/i2)(l 
(l-/.2)(l-3/.'2 3(1-/.2)(1 

Py = ^1- ll^^l- fi'^-^P^^^R 



3 
8 



(B.43) 
(B.44) 

(B.45) 
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/iyU-'C yU/i'C —fiS 




(B.46) 



(B.47) 



(B.48) 



with C = cos(0 — (J)'), S = sin(0 — cj)') and 

Ct = cos(2(0 - 0')), = sin(2(0 - 0'))- The parts P^, P^^, Pt 
of P describe the scattering of the scalar, vector and tensor con- 
tribution to S respectively. 

The functions M, Q and U depend on the wave vector k, the 
photon direction n and conformal time r. We choose for each 
fc-mode a reference system with z-axis parallel to k. For scalar 
perturbations we achieve in this way azimuthal symmetry — the 
right-hand side of the Boltzmann equation and therefore also the 
brightness A4^^^ depend only on fi = [k-n) and can be expanded 
in Legendre polynomials. The right-hand side of the Boltzmann 
equation also determines the azimuthal dependence of vector and 
tensor perturbations. One can continue with this approach, but 
the resulting equations for Q and U and especially also their 
power spectra depend explicitly on the coordinate system. There- 
fore, we prefer an approach which is inherently covariant under 
rotation. 

B. 3 Electric and magnetic polarization 

Since the Stokes parameters explicitly depend on the coordinate 
system, and Eqs. ( [b.4U|) and ( |b.4i| ) transform in a rather compli- 
cated way under rotations of the coordinate system. A method to 
characterize CMB polarization due to non-relativistic Thomson 
scattering, which is more convenient than the Stokes parameters 
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since its transformation properties are very simple, has been de- 
veloped a couple years ago [|i;jU| , |it32| , P^ , P51 , [72| . A detailed deriva- 
tion of this method goes beyond the scope of this report. Here 
we just repeat the definitions and the main results. We set 



r={M,Q + iU,Q-iU) 



;b.49) 



In terms of this three component vector the collision integral 
above can we written (in vector form) as 



C[T] = aaTne[-T+ (^j- J dQ'M' + (n • v^), 0, oj 
+^ E / dQ'p('")(n,nOr' 



:b.5o) 



From Eqs. to (|b.at3|) one can determine the scattering 

matrix for the vector T 



( 



pirn) 



2 -^2 ~V22-^2 ^2 ~V2-2-^2 ^2 

-V6Yi"'2Yi" SiYi^'iY^ Z^^Yo'^Y^ 



-V6Yy^2Yo"' Z2Y2 -2Y2 '^-^Y^ -2Y^ 



(B.51) 



2 / 



where ^^1^ = sY^W and sY{^ = , 



Yi^*{n) are the spin- weighted 



spherical harmonics [|i 1 i| , |S5| , |76| 



We now decompose the Fourier components of the temperature 
anisotropy A4 and the polarization variables Q and U as 



£ m=-2 
£ m=-2 



(B.52) 
(B.53) 



Here m = is the scalar mode, m = ±1 are the vector and 
m = iJz2 are the tensor modes. The functions are closely 
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related to the spin weighted harmonics g^l 



,GT(n) = (-0' 



\ 



47r 



The evolution equations in term of these variables can be given 
in the following compact form [)72[ 



M 



k 



TTl Tfl 



Am) Q[m) 



2^ + 3 



E^r^ - k 



77.(m) 



2m (to) 2«^f^l j-,(m) 



ne(TTa[^)"^ + V6C("^)^£,2] > 2) 



(B.54) 



B.55) 



B 



(to) 



A; 



2^^ n(m) ^m t-,(to) 2«^^+1 j^(to) 



2^-1 

(to) 



2^ + 3 



;b.56) 



-nedTflB)'"^ > 2) . 
where we have set 

5^°) = ne(JraME)°\ sf^ = nearaWb + 4A;(* - $), 
5f ) = nedTflCW, 5^ ^ = ne(TTa4V2a;6, (B.57) 

5^^^ = fiearaC^^^ + 4v/273/cS, 5^^^ = fiearaC^^^ - 8./2/3H 



and C^"') = ]^[A^2'"^ — VG-E'I"'^]. The coupling coefficients are 



Note that for scalar perturbations, m = 0, 5-polarization is not 
sourced and we have Bf^'' = 0. 

Finally, we want to connect the intensities Al^'"^ with the more 
familiar expansion of the scalar (5), vector [V) and tensor (T) 
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contributions to the brightness function in terms of Legendre 
polynomials. Usually one sets 



Here Ai^^^ only depends on ji = {n-k)/k and the n-dependence 
of M^^^ and M.^^^ can be written as 



M^^\^JL, 4>) = [a^P(m) cos (i)+M^^\i^) sin(/)] (B.58) 

M^^\li, 0) = (1 - H^) cos(20) + M^P sin(20)] , (B.59) 

where is the azimuthal angle in the plane normal to k. Each of 
the functions Aii^'^'-^^fi) is now expanded in Legendre polyno- 
mials 



(B.60) 



The coefficients cr^*^'^'^^ are then related to Aif^^ via the identities 



M 



r = (2£ + lVf 



1 



(B.61) 



(B.62) 



(£ + 2)! 
with o-j^^^^ = d^^] T ^f^x^] 



JT) . 2(2£+l) m 



(2£- l)(2£ + 3) 



(B.63) 



We do not repeat this correspondence for the Stokes parameters 
Q and U since it is rather complicated and not very useful as it 
depends on the coordinate system chosen. 
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B.4 Power spectra 



Since they are functions on a sphere, the observed CMB anisotropics 
and polarization are conveniently expanded in spherical harmon- 
ics: 6T{n,To)/TQ = Y.£m CLimyi^ (n) . The coefficients aim are ran- 
dom variables with zero mean and rotationally invariant vari- 
ances, Ci = {\ aim The correlation function of the anisotropy 
pattern then has the standard expression: 

(^(ni)f^(n2)) = ^E(2^+ l)C.P.(cos^) (B.64) 

Mo Iq / 4:7T g 



where cos 9 = ni-n2 and (• • •) denotes ensemble average. We use 
the Fourier transform normalization 

/>) = 1^1 /(x) exp{ik ■ x) A , (B.65) 



with some normalization volume V. Using statistical homogeneity 
we have 



-(ni)— (n2) j = -Jdx (^^(x, ^2) I 

^ / d'k (^^{k, n,)^^{k, n^)) . (B.66) 



Inserting our ansatz ( |B.bU|) for ^ = ^A^, and using the addition 



theorem for spherical harmonics, 

Pi{ni ■ 712) = 2|fxEm>^4(ni)Y£m(n2), we find 



To To I StT i^i'^rn,m 



X 



k'dkd^jYl^^{k)Y,m'{k){aial){k) 
= ^j:{'^i^l)Pi{n^- 712) J k'dk{aea}){k) , (B.67) 
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from which we conclude 

Stt 



Cr^'^'^ = ^Jedk{\4'\to,k)f)^ (B.68) 



where the superscript ''^^ indicates that Eq. ( [b.b8|) gives the con- 



tribution from scalar perturbations and •^•^ means that it is the 
contribution to the intensity perturbation. 

The QQ, UU, MQ, MU and QU correlators depend with the 
Stokes parameters on the particular coordinate system chosen. 
It is much more convenient to express the polarization power 
spectra in terms of the variables E and B which are independent 
of the coordinate system. Furthermore, since B is parity odd, its 
correlators with A4 and E vanish. One finds the simple general 
expression [[72 



{2i + l)2cf = ^ / edkX^'^^Y}'^^' , (B.69) 



where rim = 1 for m = and rim = 2 for m = 1, 2, accounting for 
the number of modes. Here X and Y run over A4, E and B. 
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